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ABSTRACT 


Methods  of  obtaining  the  Hugoniot  equation  of  state 
were  investigated.  Several  of  these,  employing  high  explosive 
devices,  were  used  to  obtain  Hugoniot  data  for  mineral  samples 
common  to  many  igneous  rocks.  Hugoniot  data  were  found  for 
orthoclase,  oligoclase,  labradorite  and  olivine  in  the  pressure 
range  from  50  to  300  kilobars.  Analytical  synthesis  models 
were  constructed  and  used  to  determine  the  synthesized  Hugoniot 
equations  of  state  for  granodiorite,  gabbro  and  dunite.  These 
compared  favorably  with  existing  Hugoniot  data  for  similar 
materials.  Methods  were  also  developed  and  used  to  predict, 
roughly,  Hugoniot  curves  for  other  geological  composites  for 
which  no  experimental  data  are  presently  available.  These 
materials  included  syenite,  quartzdiorite,  diorite,  olivine 
diabase,  and  diabase.  Estimates  of  temperatures  along  the 
Hugoniots  and  along  several  selected  unloading  adiabats  were 
calculated  for  several  minerals  and  igneous  rocks. 
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I  INTRODUCTION 


In  many  applications,  knowledge  of  the  dynamic  equation  of 
state,  or  Hugoniot,  of  a  material  is  required.  Various  laboratory 
and  field  experiments  have  been  devised  to  determine  the  Hugoniot 
of  a  particular  material.  These  experiments  are,  in  general, 
costly  and  time-consuming,  while  the  range  of  materials  of  inter¬ 
est  is,  in  some  cases,  extensive.  This  is  especially  true  where 
geological  materials  are  involved.  These  exist  in  nature  in 
almost  endless  variety  depending  only  on  the  relative  abundance 
of  the  constituents  of  the  material.  Although  many  geological 
materials,  which  have  been  and  will  be  of  interest,  differ  only 
slightly  in  mineral  content,  it  has  been  customary  to  investigate 
each  composite  separately.  Thus,  previous  Hugoniot  informs -ion 
on  similar  materials  has  been  utilized,  at  best,  only  in  a 
qualitative  way. 

One  aim  of  this  program  was  to  determine  the  feasibility  of 
utilizing  icnown  constituent  iiugoniot  data  to  determine  analyti¬ 
cally  the  Hugoniots  of  -composites  made  up  largely  of  these  con¬ 
stituents.  Consideration  was  given  to  one  class  of  geological 
materials,  the  igneous  rocks.  The  Hugoniots  of  several  igneous 
rocks  have  been  previously  determined.  In  addition,  a  few  of  the 
mineral  constituents  common  to  several  rocks  have  been  investi¬ 
gated. 

Before  this  synthesis  method  could  be  studied,  it  was 
necessary  to  obtain  Hugoniot  information  for  several  additional 


minerals,  Orthoclase,  Labradorite,  Oligoclase,  Clivine  and 
Biotite  were  investigated  in  this  program.  This  choice  was 
based  on  the  abundance  of  these  minerals  in  many  igneous  rocks 
and  on  their  structural  similarity  to  other  minerals  occurring 
in  rocks  for  which  Hugoniot  information  is  presently  known  or 
anticipated. 

Many  factors  other  than  the  mineral  content  of  an  igneous 
rock  might  influence  its  Hugoniot.  These  include  grain  size  and 
orientation  of  the  constituent  minerals,  water  content,  and 
porosity.  The  way  in  which  each  of  these  factors  affects  the 
response  of  the  composite  to  dynamic  loading  would  form  the  basis 
of  an  ultimate  synthesis  theory.  The  application  of  such  a  theory 
would  greatly  reduce  the  effort  involved  in  composite  Hugoniot 
determination.  In  application,  it  would  be  necessary  only  to 
measure  the  necessary  factors  and  apply  the  analytical  expression 
relating  them  to  the  Hugoniot  of  the  composite. 

In  this  investigation,  the  synthesized  Hugoniots  were  based 
on  the  measured  or  assumed  mineral  content  only.  Although,  for 
the  mineral  data  used  in  the  synthesis,  Hugoniot  data  were  avail¬ 
able  for  several  crystal  orientations,  only  an  average  Hugoniot 
was  used  for  each  mineral.  Even  with  this  simplified  model  the 
results  are  very  promising.  Unfortunately,  for  several  of  the 
composites  for  which  Hugoniot  measurements  have  been  made,  little 
is  known  about  the  mineral  content.  Future  experiments  will 
have  to  be  performed  on  composites  of  carefully  controlled  or 
measured  structure  and  petrography. 
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Within  the  frarnework  of  the  simplified  synthesis  model, 
two  methods  were  used  to  calculate  composite  Hugoniots.  The 
first  was  direct  application  of  the  Hugoniots  of  the  known  con¬ 
stituent  minerals  to  determine  the  Hugoniot  of  the  composite. 

In  the  second,  "indirect,"  method  both  known  mineral  data  and 
known  composite  data  were  used  to  increase  the  number  of  "known" 
minerals.  These  minerals  were  then  used  to  expand  composite  cal¬ 
culations  to  those  containing  relatively  large  amounts  of  minerals 
for  which  experimental  Hugoniots  are  not  yet  available.  Compos¬ 
ites  for  which  these  "implied"  Hugoniots  were  calculated  are 
syenite,  quartzdiorite,  diorice,  olivinediauase,  and  diabase. 

These  Hugoniot  data  are  considered  to  be  of  more  value  for  illus¬ 
trative  purposes  than  for  use  in  calculations  requiring  exact 
Hugoniots  for  these  materials.  This  indirect  synthesis  illus¬ 
trates  a  method  by  which  a  maximum  number  of  Hugoniots  can  be 
calculated  from  a  minimum  number  of  experimentally  determined 
Hugoniots . 

An  attempt  was  made  to  use  only  single  crystals  of  minerals 
in  the  Hugoniot  determinations.  This  was  not  possible  for  olivine. 
Although  the  igneous  rocks  are  made  up  of  polycrystals  of  the 
minerals,  these  are  not  always  in  random  orientation  and  it  seems 
desirable  to  have  single  crystal  data  for  several  orientations. 

Such  data  can  then  be  averaged,  as  was  done  here,  or  used  in  a 
more  complete  theory  where  crystal  orientation  is  considered  to 
be  a  factor. 

Many  of  the  minerals  obtained  in  single  crystal  form  were 
found  to  have  relatively  large  variations  in  density.  Because 
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of  this  variation,  it  was  decided  to  perform  the  shock  experiments 
on  only  a  small  number  of  single  crystals  of  each  type,  making 
up  all  required  samples  from  these  few.  This  procedure  often 
imposed  the  condition  of  carrying  out  the  experiments  on  very 
small  samples.  It  was  necessary  to  examine  some  of  the  experi¬ 
mental  techniques  previously  used  and  some  variations  of  these 
techniques  to  determine  the  feasibility  of  each  for  the  intended 
application.  This  examination  formed  one  of  the  three  major 
phases  of  this  investigation. 

The  other  two  phases  were  devoted  to  the  experimental  de¬ 
termination  of  the  mineral  Hugoniots  and  the  synthesis  of  com¬ 
posite  Hugoniots  from  these  results.  Discussions  of  the  results 
and  methods  used  in  each  of  these  phases  constitute  the  major 
sections  of  this  report.  The  final  section  summarizes  the  major 
results  of  the  study  and  briefly  outlines  the  future  requirements 
and  outlook  for  Hugoniot  synthesis. 
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II  METHODS  OF  HUGONIOT  DETERMINATION 

A.  Hugoniot  Description 

I.  Shock  Wave  Equations 

The  dyr.aniic  loading  process  gives  rise  to  the  dynamic  equa¬ 
tion  of  state, or  Hugoniot,  of  a  material.  In  dynamic  loading, 
diffei'ent  regions  of  the  material  are  at  different  stress  levels 
and  have  correspondingly  different  values  of  density  and  temper¬ 
ature.  For  this  reason,  although  local  thermodynamic  equilibrium 
is  supposed  to  exist,  the  whole  sample  is  not  in  equilibrium. 

Each  region  of  local  equilibrium  anu,  therefore,  uniform  tempera¬ 
ture,  pressure,  and  density  corresponds  closely,  in  the  absence 
of  shear  forces,  to  a  state  in  hydrostatic  compression.  As  a 
result,  not  only  does  the  Hugoniot  give  basic  information  re¬ 
garding  the  shock  loading  process  but,  in  addition,  relates  to 
the  fundamental  properties  of  materials. 

In  Figure  1,  the  shock- loading  process  is  compared  with  a 
similar  static  loading  process;  both  represent  cases  of  uniaxial 
strain,  i.e.,  strain  occurring  only  in  the  axial,  1,  or  loading 
direction,  all  other  strains  being  equal  to  zero.  As  is  seen  in 
the  static  case.  Figure  1(a),  the  whole  sample  is  in  a  state  of 
equilibrium.  The  stressed  element  shown  in  the  insert  is  the  same 
throughout  the  sample.  In  the  dynamic  case.  Figure  1(b),  the 
loading  device  causes  a  shock  wave  to  be  transmitted.  In  this 
case  the  sample  is  separated  into  two  different  equilibrium  re¬ 
gions.  The  region  to  the  right  of  the  shock  wave  discontinuity 
has  not  yet  been  compressed.  The  Hugoniot  relates  the  values 
of  stress  to  strain  e  or  specific  volume  V  in  the  shocked 
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(c)  Stress-strain  elements  in  shock  compression  zones 


Figure  1  Comparison  of  Static  and  Dynamic  Loading  Conditions 


state  to  those  In  the  initial  state.  It  refers  only  to  these 
end  states  and  gives  no  information  on  the  loading  path  between 
these  end  states.  If  the  stress,  and  temperature,  T,  be¬ 

hind  the  shock  wave  corresponded  to  those  of  the  hypothetical 
static  configuration,  the  elements  of  stress  and  strain  in  both 
would  be  equal. 

For  some  materials  for  which  a  thermodynamic  equation  of 
state  is  available,  e.g.,  a  perfect  gas,  the  Hugoniot  may  be 
derived.  For  solids,  there  is  no  general  equation  of  state  that 
may  be  used  in  such  a  derivation  and  the  Hugoniot  is  usually  not 
represented  in  an  analytical  form.  There  are,  however,  general 
shock  wave  equations  which  apply  to  all  materials.  These  equa¬ 
tions  include  only  the  stress,  in  the  shock  propagation  di¬ 

rection  and  the  uniaxial  strain  in  that  direction.  The  stress 
in  this  direction  is  referred  to  as  the  pressure.  For  very  high 
stresses,  which  might  cause  the  solid  to  liquify,  the  diagonal 
stress  elements  are  equal  to  the  pressure. 

Regardless  of  the  stress  level,  the  correct  interpretation 
is  that  the  pressure  P  is  considered  equal  to  the  diagonal  stress 
tensor  element  in  the  direction  of  the  one  dimensional  strain. 

In  the  experiments  described,  the  configuration  is  such  that  all 
measurements  are  made  in  the  one  dimensional  region  (Figure  1(c)). 

Conservation  of  mass  and  momentum  across  the  shock  wave  re¬ 


quire  that 

Po^s  "  Pl^^s"  ^Pl^ 

fo  +  Po'''  =  '■l 


(Conservation  of  mass)  (1) 
(Conservation  of  momentum) (2) 
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where 


Ug  *  the  shock  velocity 

Up-,  “  the  particle  velocity  behind  the 
shock  wave 

-  initial  pressure 

«  initial  density. 

The  pressure  behind  the  shock  wave, 

'■l  -  Po^s^Pl  +  fo  • 

In  most  cases  of  interest  here,  the  shock  compression  is  of  Che 
order  of  at  least  several  thousand  atmospheres  so  that  the  ambient 
pressure  or  initial  pressure  of  about  one  atmosphere  can  be  ne¬ 
glected.  Then, 

(4) 


From  the  definition  of  strain, 

v„-  V  p  .  P, 


1-^ 


where  V  is  the  specific  voltune,  and  equation  (1), 

^  1  PO  _  ^P1  /CN 

so  that  the  dynamic  stress  and  strain,  Pj^  and  can  be  deter¬ 
mined  from  Ug  and  Up,  the  variables  usually  directly  or  indirectly 
measured  in  shock  wave  experiments.  From  equations  (4)  and  (5), 
these  velocities  can  be  related  to  the  dynamic  material  properties 
occurring  behind  a  single  shock  wave: 


f  ^1 
PoH 
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Thus,  by  measuring  any  cwo  of  the  four  variables  occurring  In 
equations  (1)  and  (3^,  the  other  two  may  be  uniquely  determined. 
The  Hugoniot  generally  refers  to  any  representation  of  two  of 
these  variables,  though  most  often  it  refers  to  the  pressure - 
density  states. 

Up  to  this  point,  as  shown  in  Figure  1(b),  only  a  single 
shock  wave  was  considered.  For  various  combinations  of  material 
properties  and  shock- loading  pressures, two  or  more  shock  waves 
may  result  (Ref.  1).  In  such  cases,  an  element  of  material  is 
successively  loaded  by  several  shock  waves.  To  determine  the 
final  state  after  the  last  loading  shock  wave  the  equations  of 
mass  and  momentxjm  may  also  be  successively  applied.  Equations 
analogous  to  (4)  and  (5)  but  for  N  waves  are 


^  ■  ^N-1 


,  „  i-N  (Usi  -  Up,) 

"  ■  i^l  -  'Jp”) 


To  calculate  the  stress  and  strain  it  is  necessary 
to  measure  the  shock  velocities  of,  and  particle  velocities  be¬ 
hind,  these  N  waves.  Such  instrumentation  techniques  become 
quite  involved  and  analysis  becomes  especially  difficult  where  N 
is  large  and  samples  are  too  short  to  obtain  the  necessary  readings 
prior  to  the  occurrence  of  shock  reflections. 
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Both  dynamic  yield  points  and  material  phase  changes 
occurring  below  the  shock-loaded  high  pressure  state  may  cause 
multiple  shock  systems.  In  these  experiments  on  single  crystals, 
it  was  anticipated  that  both  might  be  experienced  and,  therefore, 
consideration  was  given  to  instrumentation  systems  capable  of 
yielding  the  required  information. 

It  is  often  necessary  to  deal  with  only  two  waves,  a  pre¬ 
cursor  having  the  velocity  of  a  longitudinal  acoustic  wave  (Ref.  1), 
and  a  plastic  wave  Ug2  whose  velocity  is  dependent  on  the  pressure 
P2.  In  this  case,  equations  (8)  and  (9)  become 


=  P,  + 


P.(U 


S2 


Upi)(Up2 


Upi) 


1  -  Up^/C 


(10) 


1  - 


(C  -  UpP(Us2  -  Up2) 
C(Us2  -  Upi) 


PI' 


(11) 


where  the  shock  velocity  is  taken  as  C  the  acoustic  velocity 
and  Pj  is  the  dynamic  yield  pressure. 

The  Hugoniot  for  a  material  which  would  be  expected  to  yield 
a  two-wave  system  is  shown  in  Figure  2(a).  For  such  a  material 
shocked  to  a  pressure  P25  the  resulting  pressure  pulse  is  shown 
in  Figure  2(b)  and  the  regions  of  various  instantaneous  pressures 
in  Figure  2(c). 

As  in  a  one-wave  system,  the  wave  and  particle  velocities 
are  found  from  equations  (10)  and  (11)  to  be  functions  of  the 
slope  and  area  under  the  P-e  Hugoniot  curve.  Equations  analogous 
to  (5)  and  (6)  for  a  two-wave  system  are: 
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vJhere  Up^  is  the  particle  velocity  behind  the  first  or  precursor 
wave.  The  precursor  wave  transmitting  the  pressure,  =■  Pyp, 
travels  at  a  higher  velocity  than  Ug2.  As  the  second  wave  ar¬ 
rives  at  a  point  in  the  material,  the  p;  ssure  at  that  point  is 
raised  from  P^  to  P2.  The  final  particle  velocity  behind  Ug2  is 
greater  than 

u  .  u  =  /(P;  •  fpU;  - 

^P2  ^P1  /  /  p 

V  o 

It  may  be  seen  that  the  velocity  of  the  second  wave  Ug2  is  equal 
to  the  precursor  for  f 2  ^  ^2  ^  ^YP’  Therefore,  for  this 

material,  only  shock  pressures  between  P^p  and  Pg  give  rise  to 
two-wave  sj's terns. 

For  such  a  two-wave  system,  it  would  be  necessary  tc  measure 
the  shock  velocities  of  botn  waves  C  and  Ug2  and  the  particle 
velocities  Upj^  and  Up2  behind  each  wave.  From  this  information, 
two  Hugoniot  states  are  determined; 

(Pf.ej^)  and  (F2>^2^ 

where  P^^  =  PqCUpj^  is  the  dynamic  yield  pressure  and  -  Up^/C 
is  the  strain  at  the  yield  point.  P2  and  £2  loafer  to  the  high 
pressure  state  behind  the  second  wave.  For  such  a  "aterlal  the 
complete  Hugoniot  would  be  found  from  several  shock  experiments 
at  different  final  pressures  P2.  In  each  experiment  the  values 
of  P^  and  would  be  expected  to  be  the  same.  Recently  it  has 
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been  pointed  out,  however,  that  stress  relaxation  effects  can 

I 

cause  an  initially  higher  apparent  yield  point  (Ref.  2),  de¬ 
caying  gradually  to  the  dynamic  yield  pressure  Pj .  When  these 
effects  are  present,  the  measured  value  of  the  precursor  wave 
velocity  iould  depend  on  the  sample  length.  Furthermore,  studies 
of  the  equation  of  state  of  snow  (Ref.,  3)  indicate  that  the  value 
of  the  dynamic  yield  point  is  not  independent  of  the  value  of  the 
final  pressure,  but  increases  with  P2.  Both  the  determination  of 
constituents'  Hugoniots  and  the  derivation  of  an  ultimate  synthesis 
model  would  be  complicated  by  these  features,  if  significant. 

2.  Analytical  Form  of  Hugoniot 

A  Hugoniot  with  a  single  discontinuity  at  P^  is  adequate  to 
describe  many  materials.  For  a  large  number  of  materials,  the 
Hugoniot  is  adequately  represented  by  a  simpler  curve  in  (P,e) 
coordinates,  one  in  which  P^  is  sufficiently  small  that  it  may  be 
neglected.  The  Hugoniot  then  has  the  appearance  of  that  in 
Figure  3(a).  These  materials  may  always  be  loaded  by  a  single 
shock  wave  (Ref.  1).  For  a  large  nxjmber  of  materials,  the  ex¬ 
perimentally  determined  variables  of  particle  velocity  and  shock 
velocity  associated  with  this  single  wave  are  linearly  related 
(Ref.  4),  as  Figure  3(b)  shows.  The  complete  Hugoniot  is  then 
described  by  only  two  states,  and  the  shape  of  the  Hugoniot  in 
the  (P,6)  plane  is  fixed  by  the  slope  and  intercept  of  the  (Ug,Up) 
curve.  The  resulting  Hugoniot  equation  is 

p  C^e 

P  =  — - 2  •  (12) 

(1-Se)2 


Low  yield  point  neglected 


(a)  Material  with  negligible  yield  pressure 


Ug  =  C  +  SUp 


(b)  Linear  (Ug,Up)  plot 


Figure  3  Material  Characterized  by  a  Single  Shock  Wave  (a)  and 
Represented  by  a  Linear  Velocity  Relationship  (b) 


The  value  of  t  le  intercept  C  may  be  approximated  by  the  bulk 
sound  speed j  yielding 


(13) 


If  the  bulk  compressibility  is  known,  as  it  often  is,  and  if 
the  linear  relationship  between  Ug  and  Up  pertains,  it  is  neces¬ 
sary  only  to  determine  the  value  of  S  to  obtain  a  useful  descrip¬ 
tion  of  the  (P, e)  states  occurring  in  shock  compression. 

There  is  at  present,  however,  no  theoretical  justification 
for  the  linear  (Ug,Up)  relationship,  even  for  those  metals  in 
which  it  has  been  found  to  exist.  Therefore,  presently,  S  may 
be  found  only  experimentally.  For  many  metals,  however,  the 
S-values  lie  within  fairly  narrow  limits.  Using  the  data  of 
reference  4  an  average  value  of  slope  of  about  1.42  is  within 
15  per  cent  of  the  experimentally  determined  value  of  S  for  over 
80  per  cent  of  the  metals  satisfied  by  the  linear  (Ug,Up)  rela¬ 
tionship.  These  relatively  narrow  limits  suggest  that  appropriate 
values  of  S  might  be  found  for  other  classes  of  materials.  Fur¬ 
thermore,  even  for  materials  not  representable  by  a  single  straight 
line  on  the  (Ug,Up)  curve,  there  are  indications  that  a  series 
of  straight  lines  would  be  applicable  (Ref.  5,6). 

Several  minerals  are  believed  to  be  in  this  category,  i.e., 
requiring  only  two  lines,  one  below  the  yield  point  (S  =  0),  and 
one  above  the  yield  point  (S  >0).  At  present,  however,  there  is 
insufficient  mineral  data  to  determine  the  value  or  existence  of 


a  suitable  average  S. 


Equation  (13)  suggests  that,  to  a  rough  approximation,  the 
Hugoniots  for  which  an  average  S  exists  depend  significantly  only 
on  the  compressibility.  Since,  in  the  present  synthesis  method, 
the  (F,i)  OL  (P,V)  states  of  the  constituents  are  required,  it 
was  desired  to  further  investigate  the  significance  of  in  des¬ 
cribing  the  known  metal  Hugoniots. 

Rather  than  applying  equation  (13)  directly  with  a  particular 
value  of  S,  another  empirical  procedure  was  adopted.  The  metal 
Hugoniot  data  of  reference  4  were  replotted  against  the  zero  pres¬ 
sure  compressibility  as  listed  in  references  7  and  8  for  several 
values  of  the  strain  for  metals  described  by  the  linear  (Ug,Up) 
relationships.  These  curves,  shown  in  Figure  4,  were  used  to 
generate  semiempirical  Hugoniot  data  by  picking  the  appropriate 
V  lue  of  the  compressibility  for  a  metal  and  simply  reading  the 
pressures  for  each  of  the  isostrain  curves  of  Figure  4.  These 
are  compared  with  the  (P,£)  states  for  several  metals  investigated 
in  reference  4  in  Figure  5. 

Since  we  are,  in  reality,  only  comparing  the  derived  (P,V) 
states  with  the  original  data  from  which  they  were  derived,  the 
fairly  good  agreement  is  not  startling.  V/hat  is  more  interesting 
is  the  importance  of  as  a  single  parameter  describing  the  (P,V) 
states  for  this  class  of  materials.  The  scatter  in  the  Figure  4 
data  indicates  that  a  single  slope  S  is  not  exact  in  descri.bing 
all  of  these  metals. 

Other  Materials 

Unsuccessful  attempts  were  made  to  use  the  curves  of  Figure  4 
to  generate  Hugoniot  data  for  several  geological  materials.  Be¬ 
cause  many  geological  composites  have  bilinear  (Ug,Up)  relationships 
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Figure  5  Comparison  of  Hugoniot  Data  (Ref.  4)  and  Hugoniot 
Based  on  Average  Log  P,  Ix)g  p  Plot  of  Figure  4 


(Ref.  5,6,9,10),  the  S- values  associated  with  the  higher  pressure 
states  reached  by  a  single  shock  wave  are  concentrated  at  values 
around  1.6  with  a  somewhat  wider  spread  than  for  the  metals 
(Ref.  10).  Interestingly,  the  value  of  S  forair,  for  which  the 
Hugoniot  is  analytically  derivable,  converges  to 

,  where  y  =  Cp/Cy  =  1.4 

for  high  pressure  states.  Similarly,  for  a  monatomic  perfect 
gas,  S  converges  to  1.33,  while,  for  a  more  complicated  perfect 
gas,  S  is  lower  yet  always  greater  than  1.  An  average  linear  fit 
using  the  analytical  curve  of  reference  11  gives  a  value  of  ~1.6 
for  water.  A  large  fraction  of  the  materials  for  which  Hugoniot 
data  are  known  have  S- values  within  the  range  of  1  to  1.67.  Al¬ 
though  the  analytically  derived  Hugoniot  for  a  perfect  gas  could 
not  be  expected  to  apply  to  other  materials,  the  implied  signifi¬ 
cance  of  Y  is  notable. 

B.  Experimental  Hugoniot  Determination 

The  experimental  work  in  this  program  was  divided  into  two 
phases.  The  first  phase  was  directed  toward  an  examination  of 
various  experimental  techniques  to  determine  those  best  suited 
for  use  with  the  samples  of  interest.  The  second  phase  was  de¬ 
voted  to  determining  the  Hugoniots  for  the  minerals,  utilizing 
the  results  of  the  earlier  experimental  work.  The  results  of  both 
phases  and  the  resulting  mineral  Hugoniot  data  are  discussed  in 
this  section. 

The  minerals  of  interest  were  some  of  the  major  constituents 
of  several  igneous  rocks.  Those  investigated  were  orthoclase. 
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oligoclase,  labradorite  and  olivine.  Certain  peculiarities  of 
the  feldspar  group  made  the  complete  separation  of  these 
phases  impossible.  These  peculiarities  are  associated  with  the 
apparent  existence  of  a  very  weak  precursor  wave.  Shock  wave 
experiments  of  an  exploratory  nature  were  required  and  performed 
during  the  entire  period  of  the  investigation.  The  resulting 
interrelationship  between  the  two  phases  makes  it  convenient  to 
discuss  both  phases  together.  The  only  exception  is  the  following 
brief  discussion  of  some  of  the  early  experiments  that  were  en¬ 
tirely  of  an  exploratory  nature. 

Before  embarking  on  a  detailed  discussion  of  the  shock  ex¬ 
periments,  it  is  necessary  to  consider  the  relationships  between 
the  measurable  variables  and  the  der^ired  information  to  be  derived 
from  these  variables.  Although  many  references  concerning  these 
relationships  are  available  in  the  literature,  it  is  desirable  to 
re-emphasize  some  of  that  work  for  certain  aspects  of  the  present 
problem. 

1.  Shock  State  Determination 

The  shock  equations  discussed  earlier  indicate  the  require¬ 
ment  to  measure  the  shock  and  particle  velocity  behind  each  wave 
in  order  to  determine  the  shock  states  or  Hugoniot.  The  instru¬ 
mentation  schemes  used  here  are  capable  of  furnishing  information 
on  free  surface  motion  only.  Electrically  conducting  pins  and 
piezoelectric  crystals  (Ref.  12,13)  have  been  used  to  indicate 
the  time-of-arrival  of  a  conducting  surface  at  a  new  location: 
Changes  in  light  reflectance  with  time  of  a  mirror  placed  at  the 
free  surface  of  the  sample  have  been  used  to  iueasura  the  free 
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surface  velocity  (Ref.  14,15,16).  It  has  been  shown  that  the 
free  surface  velocity  is  a  good  approximation  to  twice  the  par¬ 
ticle  velocity  for  metals  (Ref.  16).  Free  use  has  been  made  of 
this  free  surface  approximation  for  all  of  the  optical  records. 

Pins  and  crystals  placed  at  or  near  a  free  surface  have  been  used 
to  indicate  shock  wave  arrival  time  and,  thus,  shock  wave  velocity. 
Two  basic  types  of  experiments  were  used:  the  impedance  match 
solution  incorporating  the  electrical  instrumentation  and  the 
velocity  method  incorporating  optical  instrumentation.  In  several 
experiments  both  instrxjmentations  were  included. 

Impedance -Match  Method 

The  impedance  matching  method  depends  on  the  existence  of 
known  Hugoniot  data,  both  loading  and  unloading,  for  one  or  more 
witness  materials  (Ref.  17).  In  addition,  this  method  can  be  used 
only  for  single-wave  loading  conditions.  A  possible  sample  con¬ 
figuration  is  shown  in  Figure  6(a).  Here  a  v/itness  material  of 
the  same  material  as  the  buffer  plate  is  affixed  to,  or  machined 
as  part  of,  the  plate.  The  sample  material  is  also  affixed  to  the 
buffer  plate.  From  readings  of  the  times  of  motion  of  the  three 
free  surfaces,  i.e.,  the  shock  arrival  times  t  at  the  buffer  sam- 
pie  interface,  t  at  the  witness  material  free  surface,  and  t  at 
the  sample  free  surface,  the  transit  times  and  shock  velocities 
through  the  buffer  plate  and  sample  are  determined.  From  the 
known  Hugoniot  of  the  witness  material,  the  pressure  P  and  par¬ 
ticle  velocity  Up  are  known,  as  Figure  6(b)  shows.  For  the  sample 
material,  the  slope  on  the  (P,Up)  plot  is  known  from  the  measured 
shock  velocity  in  the  sample.  From  equation  (4), 


(a)  Experimental  (b)  Impedance  matching  solution 

configuration 


(c)  t“X  diagram 


Figure  6  Schematic  of  Experimental  Configuration  (a)  Suitable 
for  Impedance  Matching  Solution  (b)  Based  on  Con¬ 
tinuity  of  Pressure-Material  Velocity  States  (2,2’) 
Shown  in  t-x  Diagram  of  (c) 
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Since  the  normal  pressure  and  material  velocity  are  continuous 
across  the  sample  buffer  interface,  the  pressure  velocity  states 
in  regions  2  and  2'  of  Figure  6(c)  are  the  same.  The  state  (2) 
of  Figure  6(c)  represents  the  state  behind  the  reflected  wave. 

The  slope  of  equation  (14)  must  then  intersect  the  locus  of  possi¬ 
ble  reflected  states  in  the  buffer  material.  Since  the  buffer 
material  is  usually  of  higher  impedance,  p  U  ,  than  the  sample 
material,  these  possible  states  represent  the  unloading  adiabat 
from  the  initial  state  (1),  i.e.,  the  dashed  line  in  Figure  6(b). 

In  this  work,  the  unloading  adiabat  has  been  taken  as  the  mirror 
reflection  of  the  known  loading  curve  of  the  buffer  and  witness 
materia  1  (either  aluminum  or  brass) . 

One  advantage  of  this  method  is  that  the  free  surface  approx¬ 
imation  is  not  required  because  it  is  not  necessary  to  measure  the 
known  particle  velocity  in  the  sample.  The  attending  disadvantages 
are  that  the  location  of  state  (1)  from  shock  velocity  measure¬ 
ments  alone  is  difficult.  It  also  may  be  argued  that  the  assump¬ 
tion  regarding  the  unloading  path  in  the  buffer  material  is  more 
severe  than  the  assumption  of  the  validity  of  the  free  surface 
approximation  (Ref.  18).  Because  of  these  difficulties,  several 
experiments  combining  the  impedance  method  with  optical  methods 
were  run  in  addition  to  the  series  using  pins  alone. 

Velocity  Method 

In  the  velocity  method,  the  shock  and  particle  velocity 
associated  with  each  wave  must  be  known.  For  materials  compressed 
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by  a  single  wave,  only  two  velocities  are  required;  For  more 
complicated  materials  where  there  are  two  or  more  waves,  the  opti¬ 
cal  methods  of  Fowles  (Ref.  14)  and  Wackerle  (Ref.  15)  are  parti¬ 
cularly  suited  as  they  are  apable  of  furnishing  a  continuous 
history  of  the  free  surface  motion.  In  our  experiments, the  opti¬ 
cal  methods  were  restricted  to  several  variations  of  the  inclined 
mirror  method  of  Fowles  ('Ref.  14).  Some  laboratory  simulation 
experiments  were  made,  however,  using  the  wire  method  of  Wackerly 
(Ref.  15).  (These  experiments  are  discussed  in  Appendix  A.)  Al¬ 
though  a  few  experiments  we.  e  run  using  crystals  and  pins  to  meas¬ 
ure  the  free  surface  velocity,  optical  methods  were  found  more 
suitable  to  the  samples  of  interest. 

In  all  experiments  in  which  a  measurement  of  the  free  surface 
velocity  was  made,  use  of  the  free  surface  approximation  was  re¬ 
quired  to  obtain  the  necessary  particle  velocity.  A  more  subtle 
approximation  is  also  implied  in  all  Hugoniot  methods  utilizing 
free  surface  motion  measurements.  In  Figure  7,  a  schematic  of  a 
typical  inclined  mirror  experiment  is  shown.  Light  from  an  intens 
HE  light  source  is  reflected  from  the  mirrors  mounted  as  shown. 

The  reflected  light  from  the  mirrors  passes  through  the  slit  of 
a  rotating  mirror  camera*  and  is  focussed  on  the  film  plane.  The 
view  through  the  slit  (Figure  7(b))  is  recorded  on  the  film  where 
position  is  proportional  to  time.  Figure  7(c)  shows  an  idealized 
record.  As  the  shock-accelerated  free  surface  engages  the  re¬ 
flecting  surfaces  of  the  mirrors,  the  reflected  light  intensity 

■jV 

Beckman-Whitley  Model  189  with  special  telescopic  lens  and  streak 
attachments . 
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changes  (usually  reduced)  and  the  times  of  these  are  recorded. 

The  flat  mirrors,  located  on  the  buffer  plate  and  sample  free 
surfaces,  give  the  shock  time-of-arrival .  The  reflectance  as  a 
function  of  time  of  the  inclined  mirror  indicates  the  free  surface 
velocity  of  the  sample.  If  there  is  only  one  wave,  the  velocity 
angle  a  is  a  constant.  For  a  material  shocked  above  its  yield 
point  in  a  stress  region  where  it  is  compressed  by  two  waves,  a 
record  as  shown  in  Figure  7(d)  might  be  expected.  At  time  t2,  the 
pressure  reaches  the  mirror- sample  interface  and  accelerates  the 
rate  of  reflectance  change. 

This  shock  wave  configuration  is  shown  in  Figure  8.  Here  we 
consider  a  t-x  diagram  for  a  condition  yielding  an  elastic  pre¬ 
cursor  and  a  slower  plastic  wave.  In  Figure  8(a),  the  simple  case 
is  considered.  Here,  no  interactions  are  taken  into  account.  The 
free  surface  being  observed  is  simply  presumed  to  be  accelerated, 
impulsively  with  the  arrival  of  each  wave.  These  velocities  are 
related  to  the  dynamic  yield  point  Pyp  and  the  final  pressure 
by  the  equations  set  down  earlier.  Using  pir.s  in  such  an  experi¬ 
ment  would  require  two  arrays  of  pins,  one  for  each  free  surface 
velocity. 

A  more  detailed  view  of  the  situation  is  that  of  Figure  8(b). 
Here  it  is  seen  that  the  plastic  wave,  in  fact,  never  does  reach 
the  free  surface.  Instead,  succeeding  free-surface  accelerations 
are  caused  by  the  reverberating  elastic  wave.  While,  according 
to  the  theory,  the  final  free-surface  velocity  is  approximately 
twice  the  material  velocity  in  region  (2)  behind  the  plastic  wave, 
the  transition  to  the  final  free-surface  velocity  is  not  abrupt. 
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Time 


Figure  8(a)  Simple  Double  Shock  System  Involving 
No  Interactions 


Figure  8(b)  More  Detailed  Shcck  System  Showing  Shock  VJave 
Reverberations  and  Transition  Time  to  Final 
State 
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It  occurs  over  a  series  of  jumps  of  2n  times  the  particle  veloc- 
ity  in  region  (1),  where  n  is  the  n^  increment  of  free-surface 
acceleration.  The  time  period  between  successive  accelerations 
is  dependent  on  the  properties  of  the  sample  material  as  these 
affect  the  ratio  of  the  elastic  precursor  to  the  plastic  wave  ve¬ 
locity,  These  have  been  worked  out  for  several  cases  and  are 
listed  in  Table  I.  The  sum  of  these  times  is  the  total  transition 


Table  I 

TRANSITION  TIMES  ASSOCIATED  WITH  SHOCK  TRANSITIONS 


R  =  u^/Up 
(jisec) 

Pyp 

(kbars) 

6tj^ 

6t2 

H 

e=2 

0.5 

50 

0.75 

0.48 

0.48 

1.722 

1.682 

0.75 

50 

0.32 

0.116 

0.046 

0.055 

0.217 

0.5 

25 

0.71 

0.25 

0.285 

0.156 

0.70 

0.75 

25 

0.3 

0.07 

0.02 

0.011 

0.101 

time  necessary  for  the  free  surface  velocity  to  be  an  indication 
of  the  plastic  pressure.  Since  we  are  interested  only  in  the 
transition  time  for  the  initial  free  surface  velocity  to  the  final 
free  surface  velocity,  6t^/tj^  cannot  be  included  in  this  sum.  It 
is  necessary,  therefore,  in  some  cases,  to  require  that  the  instru¬ 
mentation  technique  be  capable  of  extremely  good  time  resolution 
in  order  to  observe  free  surface  velocity  changes.  In  fact,  the 
simpler  situation  depicted  in  Figure  8(a)  only  approximately  ap¬ 
plies  in  cases  where  the  time  period  6t  is  small  compared  with 


28 


t^,  the  time  required  for  transit  of  the  first  wave  across  the 
sample.  For  an  intermediate  case  where  t  is  too  large  to  neglect, 
a  time  resolution  of  at  least  0.1  psec  is  required.  A  transit 
time  t-^  of  betV7een  0.5  and  1.0  psec  is  typical  in  these  experi¬ 
ments  . 

The  details  discussed  here  have  not  yet  been  experimentally 
observed  with  sufficient  accuracy  for  complete  analysis.  An  at¬ 
tempt  to  observe  these  details  requires  very  high  camera  speeds-- 
such  speeds  must  be  consistent  with  the  ultimate  time  resolution 
prescribed  by  the  slit  width.  In  addition  high  camera  speed  re¬ 
sults  in  very  small  angles  in  Figure  7(c)  and  7(d).  These 
small  angles  make  it  very  difficult  to  determine  the  time  t2> 

Consideration  of  these  details,  therefore,  indicates  that 
the  transition  time  between  the  free  surface  vc-.lucit)-  associateu 
with  an  elastic  precursor  (and  characterized  by  the  angle  of 
Figure  7(d))  and  the  free  surface  velocity  associated  with  the 
final  plastic  state  (and  characterized  by  the  angle  a2  of  Figure 
7(d))  cannot  be  a  unique  time  t2.  Rather,  this  transition  must 
be  spread  over  a  period  of  time,  as  indicated  in  Figure  8(b) 

(i.e.,  E6t).  Then,  a  rounded  record  like  that  of  the  dashed  line 
in  Figure  7(d)  might  be  expected. 

Almost  all  records  for  the  feldspars  appeared  either  rounded 
as  in  Figure  7(d)  or  straight  as  in  Figure  7(c).  Certain  other 
peculiarities  discussed  in  the  next  section  led  to  a  decision  to 
treat  the  records  as  if  they  were  the  result  of  compression  by  a 
single  shock  wave. 

One  difficulty  in  attempting  to  analyze  these  records  is  the 
fact  that  any  observed  rounding  may  be  the  result  of  material 
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propercies.  Consider,  for  example,  a  material  characterized  by 
the  Hugoniot  of  Figure  9(a).  'Such  a  material  shocked  by  a  flat- 
topped  pressure  pulse  of  amplitude  P2  would  cause  the  transmittal 
of  a  pressure  pulse  of  the  type  shown  in  Figure  9(b).  Because  of 
the  dispersive  nature  of  this  material,  the  pressure  pulse  world 
change  shape  with  time  as  shown  in  Figure  9(b),  the  lower  ampli¬ 
tude  portions  of  the  wave  traveling  at  higher  velocities  than  the 
higher  amplitude  portions.  The  arrival  of  such  a  wave  at  a  free 
surface  would  cause  the  gradual  acceleration  of  the  free  surface 
rather  than  the  abrupt  acceleration  often  associated  with  a  sharp 
yield  point.  This  gradual  acceleration  would  cause  a  streak  rec¬ 
ord  of  an  inclined  mirror  to  have  a  rounded  appearance,  which,  in 
many  cases,  would  be  difficult  to  distinguish  from  the  rounding 
caused  by  shock  reverberations  at  the  free  surface. 

Records  of  this  type  are  discussed  in  the  next  section.  In 
general,  however,  while  such  records  may  be  analyzed  graphically, 
a  computer  solution  to  the  problem  is  considered  much  more  amena- 
ble  for  the  records  of  interest  here.  The  short  times  associated 
with  any  curvature  make  a  graphical  solution  too  crude.  Unfor¬ 
tunately,  no  such  analysis  was  possible  on  any  records  in  which 
curvature  was  noticed  or  suspected. 
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(a)  Hypothetical  Hugoniot  of  material  that 
would  cause  compression  fan 


(b)  Transmitted  pressure  pulse  associated  with  a 
compression  fan 


Figure  9  Material  Properties  Leading  to  Compression  Fan  of 
Shock  Waves  Resulting  in  Continuous  Change  in 
Free  Surface  Velocity 


Ill  SHOCK  WAVE  EXPERIMENTS 


In  the  experimental  portion  of  this  program,  60  experiments 
were  run.  About  one- third  of  these  were  primarily  exploratory 
and  not  intended  for  obtaining  Hugoniot  data.  Many,  however,  in¬ 
corporated  mineral  samples.  These  earlier  experim.ents  were  used 
to  investigate  the  effect  of  various  optical  parameters,  such  as 
optical  configuration,  slit  width,  light  intensity,  mirror  con¬ 
figuration  and  magnification.  In  addition,  in  many,  electrical 
time-of-arrival  indicators  were  also  used  and  comparisons  between 
piezoelectric  crystals  and  pins  were  made. 

As  a  result  of  these  preliminary  experiments,  several  modi¬ 
fications  to  the  existing  optical  system  were  decided  upon.  These 
involved  the  construction  of  a  new  lens  system  and  a.  modified  slit 
arrangement.  As  an  alternative  to  carrying  out  additional  experi¬ 
ments  of  a  purely  exploratory  nature,  a  laboratory  optical  mock-up 
was  set  up  and  used  to  compare  the  wire  technique  of  Wackerle 
(Ref.  15)  with  the  inclined  mirror  technique  for  various  free 
surface  motions.  The  results  of  this  investigation  indicated  that 
the  differences  in  sensitivity  of  the  two  methods  to  changes  in 
free  surface  motion  did  not  warrant  the  additional  effort  involved 
in  providing  the  reflecting  surface  necessary  for  the  wire  method. 
Consequently,  only  the  inclined  mirror  optical  method  was  used  in 
the  high  pressure  experiments. 

These  mock-up  experiments,  as  well  as  preliminary  high  pres¬ 
sure  experiments  on  the  minerals,  indicated  the  degree  of  diffi¬ 
culty  that  night  be  encountered  in  locating  any  times  associated 
with  the  arrival  of  a  second  or  plastic  wave  in  yielding  materials 
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(Figure  7).  The  results  of  a  brief  study  to  determine  the  effect 
of  arbitrarily  choosing  such  transition  times  are  discussed.  The 
major  data-gathering  high-pressure  experiments  were  divided  into 
several  series  involving  different  combinations  of  optical  and 
electrical  configurations. 

A.  Inclined  Mirror  Experiments 

This  series  of  experiments  utilized  the  mirror  configuration 
shown  in  Figure  10.  A  4- in.  plane- wave  generator  in  direct  con¬ 
tact  with  a  4-in.— diam  by  1/2-in, -thick  aluminiim  buffer  plate  was 
used  as  the  shock  wave  driver.  The  first  experiment  in  this  series, 
using  the  mirror  configuration  shown  in  Figure  10(a),  indicated 
the  presence  of  a  very  low  amplitude  precursor  wave.  There  was 
no  measurable  free  surface  velocity  associated  with  this  wave,  as 
may  be  seen  in  Figure  11,  so  that  the  pressure  level  transmitted 
by  this  wave  was  considered  negligibly  small.  This  inference  in¬ 
dicated  that  the  impedance  match  solution  discussed  might  be  appli¬ 
cable  to  the  minerals  of  Interest  here.  Consequently,  several 
additional  experiments  were  carried  out  using  this  mirror  configu¬ 
ration  but  incorporating  witness  materials  in  place  of  a  second 
mineral  sample.  The  impedance  match  solution  indicated  that  the 
hypothesis  of  the  existence  of  the  low-pressure  wave  was  correct 
but,  in  addition,  these  records  showed  that  there  was,  in  fact, 
some  very  slight  free  surface  velocity  imparted  after  the  arrival 
of  the  weak  precursor.  This  "rounding"  in  the  record  was  attri¬ 
buted  to  the  reverberations  of  the  precursor  between  the  free 
surface  of  the  sample  and  the  high  pressure  wave.  This  phenomenon 
was  considered  in  the  discussion  of  Figure  7. 
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inclined  mirrors  for 
free  surface  velocity 
measurement 


flat  mirror  for 

time-of-arrivalN^ 

measurement 

I  1 


Aluminum  buffer  plate 


Shot 

Sam 

1 

pie 
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Explosive 

14 

Oligoclase 

Orthoclase 

PWG(P-40) 

1  /■  1  ■» 

ID-  1  / 

Oligoclase 

A1 

PWG(P-40) 

18-19 

Labradorite 

A1 

PWG(P-40) 

Figure  10(a)  Mirror  Configuration  Used  for  Shock  and  Particle 
Velocity  Determination  for  Listed  Feldspar 
Experiments 
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inclined  mirrors  for 
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.measurement 


Shot 


20-21 


Aluminum  buffer  plate 


Sample 


Oligoclase 


Oligoclase 


Explosive 


PWG(P-40) 


PWG(P-40) 


*large  flat  mirrors  only 


Figure  10 (b)  Mirror  Configuration  Used  for  Shock  and  Particle 
Velocity  Determination  for  Listed  Feldspar 
Experiments 
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flat  mirror  for 
time-of-arrival 
measurement  > 


inclined  mirror  for 
free  surface  velocity 
measurement 


Aluminum  buffer  plate 


Shot 

Sample 

Explosive 

24 

Aluminum 

PWG  (P-40) 

25 

Labradorite 

PWG  (P-40) 

26 

Oligoclase 

PWG  fP-40) 

BHii 

Labradorite 

PWG +  2- in. PBX* 

28 

Oligoclase 

PWG  +  2- in.  PBX* 

29-30 

Orthoclase 

PWG +  2- in.  PBX* 

31-32 

Orthoclase 

PWG +  2 -in.  PBX* 

33-39 

(see  Table  II) 

(see  Table  II) 

■*2~in.PBX  refers  to  use  of  i  2  x  4-in.-diam  slab 
of  PBX  us^d  with  PWG  (P-40  plant'  wave  gonera-t^r) 

Figure  10(c)  Mirror  Configuration  Used  for  Shock  and  Particle 
Velocity  Detiirmination  for  Listed  Feldspar 
Experiments 
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Because  of  the  surprisingly  high  velocity  of  the  precursor 
wave  (greater  than  7inin/psec)  in  the  oligoclase  sample,  the  con¬ 
figuration  was  modified  to  that  shown  in  Figure  10(b).  It  was 
then  possible  to  measure  two  values  for  the  precursor  velocity. 

The  unexpected  result,  as  shown  in  Table  II,  is  that  rather  than 
showing  a  reduction,  the  average  value  of  the  precursor  velocity 
over  the  second  portion  of  the  oample  is  higher  than  that  over  the 
first  portion. 

In  addition,  there  was  no  indication  of  any  unloading  from 
the  explosive.  It  was  decided,  therefore,  that  the  driver  system 
using  a  plane -wave  generator  as  the  only  explosive  was  a  fairly 
good  system  for  our  purposes. 

The  next  series  of  experiments  utilized  the  mirror/sample 
configuration  of  Figure  10(c).  In  this  series,  relatively  large 
flat  mirrors  were  used  in  an  effort  to  determine  the  planarity 
and  tiroes  associated  with  the  wave  system.  A  sample  record  is 
shown  in  Figure  12.  The  results  of  these  experiments  are  also 
listed  in  Table  II.  For  several  experiments,  two  wave  velocities 
are  shown  but  only  the  particle  velocity  associated  with  the  second 
wave  has  been  listed.  The  pressure  and  strain  are  calculated  from 
equations  (4)  and  (5) . 

The  same  optical  configuration  was  used  in  an  additional 
series  of  experiments  that  incorporated  piezoelectric  crystal  trans¬ 
ducers  of  the  type  described  in  reference  12.  These  gages  are 
shown  in  place  on  a  sample  in  Figure  13.  It  was  desired  to  com¬ 
pare  the  shock  velocities  derived  from  these  gage  readings  with 
those  obtained  from  the  mirrors. 
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Figure  12  Sample  Record  of  Configuration  Shown  in  Figure  10(c)  Showing 
Arrival  Time  of  Shock  Wave  of  Flat  Mirrors  t  ,  Transit  Times 
of  Precursor  and  2nd  Shock  Waves  Atj^  and  A  to,  and  Free 
Surface  Angle  a 


Table  II 


RESULTS  OF  SHOCK  WAVE  EXPERIMENTS  ON  FELDSPARS 
(ORTHOCLASE,  OLIGOCLASE,  AND  LABRADORITE) 


Material 

Crys¬ 

tal 

Direc¬ 

tion 

Density 

(gm/cc) 

Explosive 
PWE  H  E 

Configuration 
Opt.  Elect. 

Sample 

Length 

(in.) 

14 

Orthoclase 

C 

2.57 

P-40 

•  • 

Fig 

10a 

•  • 

0.253 

14 

Oligoclase 

C 

2.65 

P-40 

-- 

Fig 

10a 

-- 

0.182 

15 

Labradorite 

C 

2.71 

P-40 

-- 

Fig 

10a 

-- 

0.265 

15 

Labradorite 

C 

2.71 

P-40 

Fig 

10a 

-- 

0.265 

16 

Oligoclase 

C 

2.65 

P-40 

Fig 

10a 

-- 

0.179 

17 

Oligoclase 

C 

2.65 

P-40 

-- 

Fig 

10a 

-- 

0.182 

18 

Labradorite 

C 

2.71 

P-40 

-- 

Fig 

10a 

-- 

0.219 

19 

Labradorite 

C 

2.71 

P-40 

Fig 

10a 

0.371 

20 

Oligoclase 

C 

2.65 

P-40 

.  - 

Fig 

10b 

-- 

0.370 

20 

Oligoclase 

C 

2.65 

P-40 

-- 

Fig 

10b 

-- 

0.186 

21 

Oligoclase 

C 

2.65 

P-40 

-- 

Fig 

10b 

-- 

0.373 

21 

Oligoclase 

C 

2.65 

P-40 

-- 

Fig 

10b 

-- 

0.191 

22 

Oligoclase 

C 

2.65 

P-40 

Fig 

10b* 

-- 

0.413 

22 

Oligoclase 

C 

2.65 

P-40 

-- 

Fig 

10b* 

-- 

0.178 

25 

Labradorite 

C 

2.71 

P-40 

-> 

Fig 

10c 

-- 

0.175 

26 

Oligoclase 

C 

2.65 

P-40 

PBX 

Fig 

10c 

-- 

0.180 

27 

Labradorite 

C 

2.71 

P-40 

PBX 

Fig 

10c 

-- 

0.182 

28 

Oligoclase 

C 

2.65 

P-40 

Fig 

10c 

-- 

0.181 

30 

Or thoclase 

B 

2,57 

P-40 

PBX 

Fig 

10c 

-- 

0.225 

32 

Orthoclase 

B 

2.57 

P-40 

Fig 

10c 

.. 

0.251 

33 

Orthoclase 

C 

2.57 

P-40 

Fig 

10c 

cryst 

0.250 

34 

Oligoclase 

C 

2.65 

P-40 

Fig 

10c 

cryst 

0.250 

35 

Orthoclase 

C 

2.57 

P-40 

PBX 

Fig 

10c 

cryst 

0.250 

36 

Oligoclase 

C 

2.65 

P-40 

PBX 

Fig 

10c 

cryst 

0.250 

37 

Labradorite 

C 

2.71 

P-40 

Fig 

10c 

cryst 

0.250 

38 

Labradorite 

C 

2.71 

P-40 

PBX 

Fig 

10c 

cryst 

0.250 

39 

Labradorite 

C 

2.71 

P-40 

Fig 

10c 

cryst 
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Table  II  (Cont.) 


Shock  Wave  Tran¬ 
sit  Times  (fisec) 
Optical  Elect 

Shock 

Velocities 
(mm/p  sec) 

Particle  Veloc¬ 
ities  (rnm/psec) 
Optical  Elect 

Pres¬ 

sure 

P 

(kilo- 

bars) 

Strain 

e 

nist 

U.qi 

Zncf 

_ 

1 

^nd 

Up9 

Est. 

Up 

ssm 

.935 

1.23 

7.02 

5.32 

102 

.141 

.585 

.873 

-  - 

7.9 

5.91 

OH 

.86 

135 

.145 

.965 

1.22 

-- 

6.99 

5.51 

O' 

- 

.60 

89 

.109 

.965 

1.22 

-  - 

6.99 

5.51 

OH 

.60 

89 

.109 

.82 

-  - 

5.55 

.70 

103 

.126 

.78 

-  - 

5.92 

.68 

107 

.115 

.87 

1.11 

•  - 

6.40 

5.05 

0 

.76 

■04 

.15 

1.29 

1.61 

-- 

7.30 

5.85 

0 

.63 

100 

.108 

1.40 

1.83 

-- 

6.7 

5.15 

0 

.78 

107 

.151 

.72 

.90 

-- 

6.56 

5.25 

0 

.78 

108 

.148 

1.31 

1.7 

-  - 

7.25 

5.60 

0 

.74 

110 

.132 

.84 

-- 

5.80 

.74 

114 

.127 

1.47 

m  m 

7.1 

-- 

-■ 

• 

.646 

7.0 

-- 

.695 

.92 

-  • 

6.4** 

4.84 

o: 

- 

.77 

iOl 

.159 

.754 

1.06 

6 . 08"^ 

4.32 

0= 

.795 

(91) 

(.184) 

.83 

-- 

5.57 

1.78 

269 

.320 

.645 

.782 

7.15 

5.89 

0- 

1.83 

286 

3.12 

.87 

1.31 

6.60 

4.36 

OH 

.77 

86 

.176 

.84 

1.17 

-  « 

7.60 

5.45 

OH 

1.50 

210 

.275 

•  • 

1.2 

5.30 

.74 

101 

.140 

.9 

1.2 

.92 

7.09 

5.30 

0+ 

.59 

.47 

83 

.111 

1.0 

1.2 

1.02 

e.2Cf^' 

".30 

0 

1.60 

1.77 

255 

.258 

.836 

1.05 

.88 

7.59 

6.05 

0 

1.78 

1.52 

286 

.294 

1.05 

1.33 

1.26 

6.05 

5.04+'- 

0 

.69 

.55 

94 

.137 

.932 

1.07 

.89 

6.82 

5.95 

0 

1.39 

224 

.234 

1.03 

1.03 

6.10 

.76 

.52 

126 

.125 

1.04 

1.3 

1.04 

6.07 

4.88 

0= 

E 

.62 

.57 

82 

.127 

1.02 

1.34 

6.24 

4.75^ 

O' 

F 

.77 

127 

.123 

5.96 

1.88 

296 

.32 

6.32 

oh: 

.74 
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.117 

%> 

5.02 

1,57 
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.314 

1.14 

1.00 

7.52 

5.58, 

.86 

.63 

135 

.145 

.964 

1.02 

9.91 

6.20+ 

0+ 

1.08 

.93 

177 

.174 

1.27 

1  11 

5.65 

.97 

.78 

145 

J 

.172 
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Experimental  Configuration  Showing  Electrical 
and  Optical  Instrumentation 

The  sample  is  from  a  single  crystal  of  ortho- 
clase.  Microcline  particles  are  also  seen. 


As  Table  II  shows,  these  velocities  do  not  generally  coin¬ 
cide  but  the  crystals  seem  to  indicate  the  arrival  of  the  plastic 
wave  and  not  the  precursor.  A  record  showing  the  times  associated 
with  both  readings  is  shown  in  Figure  14.  The  times  are  normal¬ 
ized  to  the  arrival  of  the  shock  wave  at  the  free  surface  of  the 
buffer  plate. 

It  was  hoped  that  the  crystal  gage  readings  might  indicate 
the  arrival  of  both  the  elastic  and  plastic  waves.  As  may  be  seen 
from  the  sample  record  of  the  crystals  (Figure  15),  the  gage  be¬ 
comes  inoperative  immediately  after  the  first  indication.  This 
lends  support  to  the  view  that  the  precursor  wave  is  quite  weak. 
Although  it  might  be  concluded  that  a  higher  amplitude  elastic 
wave  might  also  break  the  crystal  gage,  the  closer  correlation  of 
crystal  signal  to  the  plastic  wave  is  contrary  to  such  a  conclu¬ 
sion. 

In  an  effort  to  further  investigate  the  amplitude  and  velocity 
of  the  precursor  v/ave,  a  series  of  experiments  utilizing  the  opti¬ 
cal  configuration  shown  in  Figure  16(a)  was  run.  A  slight  modifi¬ 
cation  is  shown  in  Figure  16(b).  The  elastic  wave  velocities 
der^-ved  from  these  records  show  considerable  variation  for  a  given 
sample,  but  appear  to  be  constant  throughout  the  width  of  the  sam¬ 
ple.  For  both  oligoclase  and  labradorite,  the  high-pressure 
shots  showed  a  lower  precursor  velocity  than  the  low  pressure 
shots.  As  is  noted  in  Figures  10  and  16,  the  higher  pressures 
were  obtained  by  adding  a  2-in.  slab  of  PBX  between  the  plane 
wave  generator  and  the  buffer  plate.  A  record  from  one  of  these. 
Figure  17,  shows  the  constancy  of  the  velocity  of  the  first  wave 
through  the  sample.  This  record  is  especially  interesting  because 
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Figure  Streak  Record  (a)  from  Mirror  Configuration  of  Figure  10(c)  for  Sample  of 

Orthoclase,  Showing  Shock  Transit  Times  Measured  Optically,  Atppt,  and 
Electrically,  Atg^,  and  Free  Surface  Velocity  Angle,  a;  and  (b;^Model 
before  Firing  As  Viewed  in  Static  Frame  of  Rotating  Mirror  Camera 


Figure  15  Oscilloscope  Trace  of  Crystal  Readings  of 
Time-of-Arrival  of  Shock  Wave  at  Buffer 
Plate  Free  Surface,  t^^,  and  Sample  Free 
Surface,  t£ 
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Figure  16(a)  Mirror  Configuration  Used  for  Shock  and  Particle 
Velocity  Determination  for  Listed  Feldspar 
Experiments 
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Figure  16(b) 


Mirror  Configuration  Used  for  Shock  and  Particle 
Velocity  Determination  for  Listed  Feldspar 
Experiments 
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Figure  17  Sample  Record  oi  Configuration  16(a)  Showing  Precursor 

Shock  Velocity  Angle  aj^.  Free  Surface  Velocity  Angle  02, 
and  Lag  Time  At^  between  Arrival  of  Precursor  and  the 
Motion  of  the  Free  Surface 


it  shows  the  time  lag  between  the  arrival  of  the  first  wave  and 
the  acceleration  of  the  free  surface  wave.  Furthermore,  it  is 
seen  that  the  acceleration  of  the  free  surface  is  not  abrupt  but 
is  quite  gradual.  This  indicates  that  the  motion  is  due  to  a 
series  of  waves  generated  eithej.  from  shock  reverberations  between 
the  free  surface  and  a  later  plastic  wave  or  as  a  result  of  a 
Hugoniot  curve  rounded  in  the  vicinity  of  its  dynamic  yield  point. 

B.  Pin  Experiments 

All  experiments  described  thus  far  utilized  mirror  reflection 
techniques.  The  results  of  these  experiments  have  shown  large 
sample-to-sample  anomalies.  Nevertheless,  these  experiments  have 
been  valuable  in  gaining  insight  into  the  wave  structure  generated 
in  these  feldspars.  As  a  result  of  these  experiments,  it  was 
tentatively  concluded  that  the  high-pressure  Hugoniot  states  may 
be  treated  as  if  they  are  attained  by  shock  compression  from  a 
single  shock  wave. 

As  a  result  of  this  conclusion,  two  additional  sets  of  ex¬ 
periments  were  run.  In  the  first  set,  the  mirror  configuration 
shown  in  Figure  18  was  combined  with  the  pin  closure  switches 
shown  in  Figure  19.  The  results  of  this  series  are  shown  in 
Table  II.  The  shock  arrival  times  coincide  roughly  with  the  mo¬ 
tion  of  the  sample's  free  surface,  as  shown  in  Figure  20. 

The  second  set  of  experiments  was  based  entirely  on  the 
impedance  matching  method  and  utilized  electrical  closure  pins 
only.  In  this  set  of  experiments,  6-in.  plane  wave  generators 
were  used  with  various  explosives  and  driver  plates  to  impart  a 
series  of  loading  conditions  to  the  samples.  A  sample 
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shown  in  Figure  18  was  combined  with  the  pin  closure  switches 
shown  in  Figure  19.  The  results  of  this  series  are  shown  in 
Table  II.  The  shock  arrival  times  coincide  roughly  with  the  mo¬ 
tion  of  the  sample's  free  surface,  as  shown  in  Figure  20. 

The  second  set  of  experiments  was  based  entirely  on  the 
impedance  matching  method  and  utilized  electrical  closure  pins 
only.  In  this  set  of  experiments,  6-in.  plane  wave  generators 
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49  Oligoclase-C  dir/Alviininum  PWG(P-60)* 

50  Ollgoclase-C  dir/Oligoclase-Pdlr  PWG(P-60)* 


*6-in.-diain  plane  wave  generator 


Figure  18  Mirror  Configuration  Used  for  Shock  and  Particle 
Velocit5>^  Determination  for  Listed  Feldspar 
Experiments 
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configuration  is  shown  in  Figure  21.  In  each  experimeat,  five 
mineral  samples  were  mounted  in  five  pie-shaped  segments  of  equal 
size.  In  addition,  in  each  segment  two  witness  materials,  brass 
and  aluminum,  were  mounted.  Pin  switches  set  just  off  the  free 
surface  were  set  to  close  at  uhe  time  of  free  surface  motion. 

Pin  closure  caused  a  condenser  to  discharge  through  a  resistor 
and  the  associated  voltage  signal  was  displayed  on  an  oscillo¬ 
scope.  On  each  oscilloscope  trace,  there  were  four  signal  inputs; 
one  each  from  pins  located  to  measure  free  surface  motions  of  the 
mineral  sample,  the  two  witness  materials  and  the  buffer  plate. 

The  pin  locations  are  shown  in  Figure  22;  an  oscilloscope  trace 
is  shown  in  Figure  23.  The  pertinent  times  are  indicated  by  the 
sharp  change  in  amplitude  of  the  signal.  Fairly  good  time  resolu¬ 
tion  is  made  possible  by  a  sweep  speed  of  0.5  cm  per  psec  and  the 
signal  generator  trace  located  above  and  below  each  trace. 

The  various  minerals  of  interest  and  several  samples  of 
granite  were  distributed  among  the  eight  experiments  of  this 
series  so  that  each  sample  was  subjected  to  a  variety  of  stress 
levels.  In  Table  III  the  samples  included  on  each  shot  are  listed 
as  is  the  explosive  and  buffer  configuration. 

All  of  the  minerals  of  interest  were  Included  in  this  series 
of  experi:ments .  Samples  of  orthoclase,  labrador ite,  and  oligo- 
clase  were  cut  and  shock-loaded  in  two  directions  (referred  to 
as  B  and  C  directions).  These  orientations  are  described  in  the 
next  section.  Single-crystal  samples  of  olivine  were  not  avail¬ 
able  for  experiments  so  it  X'/as  necessary  to  use  polycrystalline 
or  "isotropic”  samples.  A  few  samples  of  biotite  were  also  in¬ 
cluded  in  this  series.  We  were  able  to  cut  biotite  only  in  one 
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(a)  Open  View 


(b)  Closed  View 


Figure  21  Configurations  Snowing  Mineral  Samples  (Square), 
Witness  Materials  (Round)  in  Open  View  (a)  and 
Pins  with  Connections  in  Closed  Vievz  (b) 
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(b)  Schematic  of  pin  locations 

Pq  and  Pq'  trigger  locations 
Pj^  buffer  free  surface 
P2-P^  0.092  in.  Al,  0.125  in.  mineral, 
0.155  in.  brass  free  surface 

Figure  22  Sample  Configuration  for  Pin  Experiments  Showing 
Mineral  and  Witness  (M^  and  Material  Loca¬ 
tions  (a)  and  Pin  Locations  P^  for  One  Segment 
and  One  Oscilloscope  Trace 


Figure  23  Oscilloscope  Trace  from  Piii  Instrumentation 
of  Figure  22,  Showing  Times  of  Arrival  at 
Buffer  Plate  Free  Surface,  t  ,  at  Aluminum 
Free  Surface,  t^;  at  Mi'icral^Free  Surface  t2 
and  at  Brass  Free  Surface,  t^  (Sweep  speed 
0.2  g,£ec/cm.  Time  marks  are'^at  intervals  of 
0.]  gsec) 


56 


Table  III 

SAMPIES  INCLUDED  IN  PIN  EXPERIMENTS 


Shot 

No. 

Exp lo  s ive 

Buffer 

Piste 

Material 

(0.500  in.) 

Mineral  Included  in 

Shot 

Orthoclase-C 

Orthoclase-B 

u 

1 

QJ 

m 

ri 

i-i 

O 

0 

bO 

•r-f 

o 

Oligoclase-P 

Labrador ite-C 

1 

o 

4J 

•H 

0 

T3 

to 

1— 1 

Olivine 

Biotite 

Granite 

_ 1 

1 

PWG 

A1 

X 

X 

X 

X 

X 

2 

P-60 

Brass 

X 

X 

X 

X 

X 

3 

P-60 

A1 

X 

X 

X 

X 

X 

4 

P-60 

Brass 

X 

X 

'X 

X 

X 

5 

P-60 

Brass 

X 

X 

X 

X 

X 

6 

P-60 

Al* 

X 

X 

X 

X 

X 

7 

P-60 

A1 

X 

X 

X 

X 

X 

3 

P-60 

Brass 

X 

X 

X 

X 

X 

*0. 125-in.  A1  flyer  plate  used;  no  data 
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direction.  Thus,  including  several  samples  of  granite,  this 
series  included  nine  sample  types.  The  results  are  listed  in 
Table  IV. 

C .  Summary  of  Experimental  Results 

The  Hugoniot  data  for  the  feldspars  are  subject  to  some  in¬ 
terpretation.  The  elastic  precursor  of  low  amplitude  is  followed 
by  plastic  loading  to  the  final  state.  The  optical  recordfi  indi¬ 
cate  that,  in  some  cases,  the  plastic  loading  is  accomplished  by 
a  single  shock  wave  while,  in  others,  the  plastic  loading  is  accom¬ 
plished  by  a  compression  fan.  Schematics  of  these  optical  records 
are  shown  in  Figure  24.  In  Figure  24(a)  the  elastic  wave  would 
be  of  very  low  amplitude  and  could  be  neglected  with  little  error. 
The  sample  is  essentially  loaded  by  the  plastic  wave  arriving  at 
t2>  In  Figure  24(b)  there  is  evidence  of  a  second  elastic  wave 
at  t'  and  a  plastic  wave  at  t2. 

The  location  of  t'  and  t2  is  very  difficult  and  these  records 
are  in  some  cases  indistinguishable  from  Figure  24(a)-  and  (c). 

In  Figure  24(c),  there  is  evidence  of  shock  loading  by  a  com¬ 
pression  fan.  Because  of  the  apparent  zero  free -surface  velocity 
in  Figure  24(c),  records  of  this  type  were  analyzed  on  the  basis 
of  a  single  loading  wave.  ‘ nterpretations  of  these  records  on 
the  basis  of  a  loading  fan  would  require  the  type  of  analysis  in¬ 
dicated  in  Appendix  B,  where  a  gradually  increasing  loading  his¬ 
tory  is  treated  as  if  a  series  of  plastic  waves  causes  the  loading. 
Although  such  an  interpretation  of  all  of  the  optical  records. 
Figure  24(a),  (b)  and  (c),  might  cause  considerable  changes  in 
the  final  results,  no  such  analysie  has  yet  been  attempted.  Each 
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of  the  records  was  considered  as  if  one  plastic  wave  arriving  at 
t2  resulted  in  the  final  state  manifested  in  the  final  free  sur¬ 
face  velocity.  This  interpretation  was  consistent  with  the  im¬ 
pedance  match  (the  pin)  experiments  valid  only  for  single  wave 
sample  loading. 

At  first  it  was  anticipated  that  the  double  shock  wave  sys¬ 
tem  ordinarily  associated  with  a  Hugoniot  with  a  sharp  yield 
point  would  lead  to  a  sharp  discontinuity  in  the  observed  free- 
surface  velocity.  As  mentioned  in  the  discussion  of  Table  1,  the 
sharpness  of  the  transition  is  a  function  of  the  separation  be¬ 
tween  the  two  waves.  In  any  case,  some  rounding,  associated  with 
the  elastic  wave  reverberations,  would  be  expected.  This  round¬ 
ing  was  especially  apparent  in  the  record  of  Figure  17 .  For  the 
same  conf i.guration,  a  slower  camera  speed  would  tend  to  give  the 
impression  of  a  sharp  discontinuity.  On  the  other  hand,  a  higher 
camera  speed  combined  with  greater  light  Intensity  and  perhaps 
a  narrower  slit  width  might  enhance  the  time  resolution  to  the 
point  where  the  reverberation  details  are  apparent.  Neither  ex¬ 
treme  was  obtained  in  this  study. 

The  methods  of  analysis  adopted  for  the  reduction  of  the 
optical  records  are  admittedly  a  practical  compromise  between  the 
rather  arbitrary  choice  of  an  intermediate  arrival  time  for  the 
plastic  wave  and  the  more  general  graphical  iteration  shown  in 
Appendix  B.  The  choice  of  the  arrival  time  of  the  plastic  wave 
between  the  two  extremes,  the  arrival  time  of  the  elastic  wave 
and  the  time  at  which  the  final  free  surface  velocity  is  achieved, 
has  profound  influence  on  the  calculated  yield  point.  The 
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graphical  method  was  attempted  on  several  records  where  rounding 
was  noticed.  Even  the  graphical  method  was  unsatisfactory  for 
the  times  and  free  surface  velocities  involved  here.  A  simple 
computer  calculation  would  be  best  but  was  not  attempted  on  these 
records . 

1.  Hugoniot  Gurves 

The  results  of  the  optical  and  combination  optical  and  elec¬ 
trical  experiments  (Table  II)  and  of  the  pin  experiments  (Table  IV) 
are  combined  into  the  Hugoniot  curves  plotted  in  pressure-strain 
coordinates  in  Figures  25,  26,  27  and  28.  In  these  figures  data 
are  listed  for  two  crystal  orientations.  These  orientations  are 
the  same  as  those  listed  in  Figures  10  and  16  and  Tables  II,  III, 
and  IV. 

In  the  Hugoniot  curves  for  orthoclase,  labradorite  and  oli¬ 
vine,  a  yield  point  is  depicted.  The  location  of  this  yield 
point  has  been  inferred  from  only  one  record  thus  far.  The  method 
of  analysis  described  obviously  precludes  the  possibility  of  such 
yield- point  determination.  Nevertheless,  Hugoniot  data  for  com¬ 
posites  made  up  of  these  minerals,  notably  granite  (Ref. 9,  10,  12), 
indicate  the  presence  of  a  yield  point  at  about  50  kilobars.  The 
yield  points  for  orthoclase,  labradorite  and  olivine  are  based 
on  a  rough  approximation  using  the  hydrostatic  compressibility 
data  of  Bridgman  (Ref.  19)  and  the  low  pressure  shear  data  of 
reference  20. 

All  of  the  synthesis  work  was  based  on  the  experimental 
Hugoniots  of  Figures  25,  26,  27  and  28.  It  is  hoped  that  a  more 
complete  future  analysis  of  the  data  could  be  based  on  a 
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Figure  25  Hugoniot  Curve  for  Orthoclase  Showing  Average  Used  in  Synthesis  Calculations 


300 
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Figure  26  .  Hugoniot  Curve  for  Labradorite  Showing  Average  Used  in  Synthesis  Calculations 


Figure  27  Hugoniot  Curve  for  Oligoclase  Showing  Average  Used  in  Synthesis  Calculations 


Normal  pressure,  P  (kilobars) 


computerized  solxition  of  the  multiwave  structure  depicted  by 
Figure  24.  It  is  noted  graphically  that  the  Hugoniots  in  figures 
show  yield  at  about  50  kilobars.  This  interpretation  is  the  re¬ 
sult  of  synthesis  calculations  on  measured  bulk  compressibilities 
and  shear  moduli  but  not  directly  on  the  measured  data. 

2.  Crystal  Orientation 

In  Figure  29,  a  crystal  of  orthoclase  is  shown.  Although 
orthoclase  is  monoclinic  while  the  plagicclases  are  triclinic,  the 
two  are  similar  because  the  triclinic  angles  between  the  a  and  b 
and  b  and  c  crysual  axes  are  close  to  90  deg.  Therefore,  in  de¬ 
fining  orientation.  Figure  28  may  be  considered  applicable  to  all 
of  the  feldspars  studied.  All  have  good  cleavage  parallel  to  the 
c-pinicoid  which  lies  in  the  plane  of  the  a  and  b  crystal  axes. 
This  cleavage  plane  was  easily  recognized  in  all  of  our  samples. 

In  addition,  many  of  our  orthoclase  samples  were  cut  from  a  sin¬ 
gle  crystal  of  orthoclase  in  which  all  of  the  pinicoids  and  crys¬ 
tal  directions  were  fairly  easily  determined.  ConseTquently,  all 
feldspar  samples  included  th^'se  cut  parallel  to  the  c-pinicoid. 
These  samples  were  mounted  so  that  the  shock  propagation  direc¬ 
tion  was  perpendicular  to  the  c-pinicoid,  as  shown  on  the  left  in 
Figure  29.  All  samples  oriented  in  this  way  are  listed  in  the 
data  tables  as  the  C-dlrection. 

For  orthoclase,  samples  were  also  cut  parallel  to  the 
b-pinicoid  (Figure  29):  These  are  listed  as  B-orientations .  Sam¬ 
ples  of  labradorite  and  oligoclase  were  not  of  a  size  or  crystal 
completeness  that  facilitated  identification  of  any  but  the 
c-pinicoids.  These  numerals  were  cut  and  mounted,  therefore,  so 
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c~pinicoid  -i 


angle  >  90 


lKp(o(  arbitrary 
but  l^p  in 
ab  plane) 

(b)  Crystal  axes 


b-pinicoi.d 


(a)  Monoclinic  crystal 


Figure  29  Monoclinic  Crystal  (a)  Slicwing  b-  and  c-pinicoids 
and  Monoclinic  Axes  (b)  Showing  Shock  Propagation 
Directions  1?!r,  and  l^p  for  Crystal  Orienta¬ 
tions  Listed  in  Tables  II  and  III  as  C,  B,  and 
P  Respectively 
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that  the  shock  propagation  direction  was  in  the  plane  formed  by 
the  a  and  b  crystal  axes,  but  otherwise  arbitrary.  Such  orienta¬ 
tions  are  labeled  with  the  symbol  P  to  indicate  the  parallel 
direction.  This  arbitrariness  may  indeed  be  partially  responsi¬ 
ble  for  some  of  the  scatter  in  the  observed  values  of  the  elastic 
wave  ^'ilocities  in  these  samples.  Ultrasonic  measurements  have 
shown  that  in  crystals  of  labradcrite  ana  oligoclase  the  values 
of  the  longitudinal  elastic  velocities  are  significantly  dependent 
on  the  crystal  orientation  (Ref.  21).  Those  values  are  not  repro¬ 
duced  here  because  of  difficulties  in  determining  the  meaning  of 
the  orientations  used  in  that  work. 

3.  Yield  Point  Approximation 

The  preceding  discussion  of  the  Hugoniot  data  and  the  inter¬ 
pretation  of  that  data  are  not  consistent  with  the  location  of 
the  yield  points  shown  in  Figures  25,  26  and  27.  Nevertheless, 
certain  other  factors  must  also  be  considered.  The  basic  objec¬ 
tive  in  obtaining  these  mineral  data  is  to  test  the  synthesis 
hypothesis.  The  feldspars  and  olivine  chosen  foi  investigation 
are  all  constituents  of  granite,  which  is  believed  to  exhibit  a 
dynamic  yield  point  at  about  50  kilobars  (Figure  30).  It  is 
reasoned  that  this  yield  point  roust  be  manifested  as  a  result  of 
similar  yielding  of  its  mineral  constituents,  individually  or  in 
combination  in  the  granite  composite.  Furthermore,  the  locatior 
of  these  yield  points,  which  is  somehow  not  apparent  from  the 
shock  data,  is  implied  from  extrapolation  of  the  hydrostatic  com¬ 
pression  data  of  Bridgman  (Ref.  19)  and  ultrasonic  measurements 
of  the  longitudinal  velocities  of  several  of  the  crystals  (Ref.  21) . 
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The  longitudinal  velocity  in  a  direction  perpendicular  to 
the  c-pinicoid  was  found  to  be  7.53  mm  per  psec  for  labradorite, 
and  6.8  mm  per  psec  for  oligoclase  (Ref.  21).  In  a  direction 
perpendicular  to  the  normal  to  the  c-pinicoid,  the  average  veloc¬ 
ity  was  6.90  mm  per  psec  for  labradorite  and  6.77  mm  per  psec  for 
oligoclase.  These  velocities  are  associated  with  linear  elastic 
stress-strain  lines  in  Figures  26  and  27.  These  lines  are  of 
slope  P/e  =  PqC  ,  where  C  is  the  appropriate  elastic  velocity. 
These  lines  of  constant  slope  have  then  been  extended  to  the  ex¬ 
trapolated  lower  end  of  the  experimental  Hugoniot  curve. 

The  hydrostatic  compression  data  of  Bridgman  for  orthoclase, 
labradorite  and  olivine  have  been  plotted  in  Figures  25,  26  and 
28.  To  roughly  establish  a  yield  point  for  these  materials,  these 
hydrostatic  data  were  used  as  a  basis  for  calculating  the  equiv¬ 
alent  normal  stress  associated  with  shock  or  uniaxial  strain.  For 
further  convenience  this  calculation  was  based  on  the  isotropic 
stress-strain  relationship  (Ref.  1) 


Tfj  »  2pe^j  +  +  £22  +  ^33^^ij 


where  t. .  is  the  stress  tensor  and  e^.  the  strain  tensor.  Al- 
though  the  LariTe  constants  X  and  p  can  be  associated  with  isotropic 
materials  only,  there  is  an  equivalent  X  and  p  for  a  polycrystal¬ 
line  material  of  lower  symmetry.  As  in  reference  1,  equation  (15) 
may  be  used  to  establish  the  relationship  between  the  hydrostatic 
and  uniaxial  compression  states: 


^l-d  =  ^hvd  +  ^ 
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For  values  of  fi.  taken  from  Birch  (Ref.  8)  this  correction  has 
been  applied  to  the  hydrostatic  data.  The  yield  points  estimated 
in  this  way  are  intended  to  give  an  average  value  of  the  Hugoniots 
of  the  minerals,  and  apply  only  to  a  quasi-isotropic  material  made 
up  of  many  crystals  in  intimate  contact  and  random  orientation. 

Each  of  the  mineral  Hugoniots  in  Figures  25,  26,  27  and  28 
consists  of  two  parts, 

(1)  an  upper  (dashed)  portion  representing  the  av¬ 
erage  between  the  two  orientations  investigated 
experimentally,  and 

(2)  a  lower  portion  representing  the  elastic  zone 
estimated  from  either  hydrostatic  or  ultrasonic 
measurements . 

The  Hugoniots  made  up  of  these  two  parts  have  been  tabulated  and 
used  for  the  composite  synthesis  calculations. 

Obviously,  the  question  of  yielding  in  these  minerals  is  by 
no  means  settled.  Additional  experiments  in  the  region  close  to 
the  yield  point,  where  the  greatest  difference  between  the  plas¬ 
tic  and  elastic  velocities  might  be  found,  are  required.  In 
addition,  more  detailed  analysis  of  some  of  the  optical  records 
exhibiting  "rounding"  will  be  informative  and  especially  so  when 
more  low  pressure  data  are  available. 

One  obvious  question  regarding  the  interpretation  of  the 
e.wperimental  records  relates  to  the  absence  of  the  high  pressures 
associated  with  the  yield  points.  From  optical  records,  such  as 
those  of  Figure  20,  it  is  implied  that  the  upper  shock  state  is 
connected  by  a  single  chord  to  a  low-pressure  state.  Such  a 
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situation  is  shown  schematically  in  Figure  31.  If  such  a  re¬ 
presentation  pertained,  one  would  expect  the  lower  plastic  veloc¬ 
ity  to  be  associated  with  the  higher  pressure  state  as  is  implied 
by  shots  40  and  41  on  Table  II. 
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Pyp  “  Hypothetical  yield  point 
implied  from  hydrostatic 
and  ultrasonic  experiments 


Figure  31  Schematic  Representation  of  Shock  Compression 
to  State  (PoEo)*  Where  Pn  >  Pyp*  Indicating 
Dynamic  Yield  Pressure  Is  Very  Low 
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IV  SYNTHESIS  OF  COMPOSITE  HUGONIOT 


A .  Advantages  of  the  Synthesis  Method 

The  usual  methods  for  determining  the  shock- loading  or 
Hugoniot  characteristics  of  a  material  are  quite  extensive  and 
costly.  These  methods  involve  laboratory  and  fielu  experiments 
of  the  type  described  in  the  last  section.  Customarily,  each 
material  for  which  the  dynamic  properties  are  desired  is  analyzed. 
In  some  cases,  although  a  new  material  is  similar  in  composition 
to  another  material  previously  investigated,  no  real  advantage  is 
taken  of  this  similarity.  The  over-all  goal  of  the  synthesis 
program  is  to  devise  a  method  by  which  the  similarities  of  nat¬ 
urally  occurring  geological  samples  may  be  used  to  simplify  the 
job  of  determining  the  dynamic  properties. 

As  an  example,  consider  the  similarity  in  composition  of  the 
principal  types  of  igneous  rocks.  Table  V  lists  the  mineral  con¬ 
stituents  (Ref.  24) .  Each  of  the  rocks  is  made  up  of  various 
percentages  of  the  14  minerals  listed.  Ideally,  the  'synthesis 
method  would  be  applied  to  this  system  as  follows; 

(1)  Analyze  a  new  rock  to  determine  the  mass  per¬ 
centages  of  the  various  minerals  present. 

(2)  Use  the  known  properties  of  the  minerals--their 
Hugcnints--to  synthesize  the  Hugoniot  of  the 
composite.  The  advantage  of  the  synthesis  scheme 
is  that  it  would  no  longer  be  necessary  to  carry 
out  the  extensive  experimentation  for  every  new 
material . 
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APPROXIMATE  MINERAL  COMPOSITIONS  OF  PRINCIPAL  TYPES  OF  IGNEOUS  ROCKS  (Ref.  24) 
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A  major  effort  on  this  program  was  to  generate  sufficient 
mineral  Hugoniot  data  to  test  the  feasibility  of  the  proposed 
synthesis  method.  Contrary  to  the  earlier  synthesis  work  of 
Porzel  (Ref.  22)  and  Chabai  (Ref.  23),  it  was  intended  here  to 
use  known  dynamic  properties  of  the  mineral  constituents  to  pre¬ 
dict  the  properties  of  composites  (here,  igneous  rocks)  made  up 
of  these  constituent  minerals.  The  experimentally  determined 
Hugoniots  to  be  used  in  the  synthesis  are  represented  by  the 
dashed  lines  in  Figures  25-28.  The  Hugoniot  cf  a  fifth  mineral 
studied,  blotite,  was  incompletely  determined  and  not  usefi;i  in 
these  synthesis  calculations  because  of  the  low  stress  levels 
attained.  The  stress  level  or  pressure  attained  in  a  given  ex¬ 
periment  is  a  function  of  the  material  itself.  The  properties  of 
biotite  are  such  that,  for  the  loading  method  used  in  these  exper¬ 
iments,  only  ~100  kilobars  was  reached,  compared  with  pressures 
of  ~300  kilobars  reached  in  the  other  materials. 

To  test  a  synthesis  method  on  rocks  of  known  properties,  the 
Hugoniot  of  quartz  was  also  needed.  The  data  for  quartz  from 
references  15  and  25  are  plotted  in  Figure  32,  Data  were  reported 
for  three  orientations  of  quartz  and  also  for  fused  quartz.  Only 
the  average  Hugoniot  is  shown  in  Figure  32.  In  this  first  at¬ 
tempt  at  synthesis,  only  one  curve  for  each  mineral  was  required; 
the  dashed  line  in  Figures  25,  26,  27,  28  and  32  represents  the 
"average”  Hugoniot  for  the  various  orientations. 

B.  Synthesis  Methods 

Two  methods  of  synthesis  were  attempted  in  this  analysis. 

The  first,  or  "direct,"  method  applies  the  known  mineral  Hugoniot 
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Hugoniot  Cxirve  for  Quartz  Showing  Average  Used  in  Synthesis 


data  to  the  determination  of  unknown  composite  Hugoniot  deter¬ 
mination.  The  equations  used  are  those  reflecting  the  law  of 
partial  volximes.  In  the  synthesis  relationship,  the  partial 
volum.e  of  a  mineral  is  taken  to  be  the  final,  volume  after  shock 
loading  to  the  high  pressvire  state.  A  more  detailed  "^eview  of 
some  of  the  differences  between  the  static  and  dynamic  multi- 
component  model  are  outlined  in  Appendix  B.  In  the  second,  or 
"indirect,"  method,  both  known  mineral  and  composite  Hugoniot 
data  are  used  to  generate  the  Hugoniots  of  unknown  composite 
materials.  The  same  synthesis  equations  are  used  in  this  method. 


1.  Application  of  Law  of  Partial  Volumes 

For  a  series  of  N  mineral  constituents  making  up  an  igneous 
composite  of  volume  V^, 

N 


where  V  is  the  volume  of  the  n"  mineral  constituent, 
n 

Since 


M  V 

c 


where  M  is  the  composite  mass  and  v^  the  specific  volume  of  the 
composite,  and 


where  M  is  the  mass  of  the  n*'^  mineral  in  the  composite  and  v 
n  n 

its  specific  volume. 


v_  = 


N 

V 

n 


M 


n 


TT 


v„ 


n 


N 

=  Z 
n 


a  \ 
n  n 


(17) 
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where  =  M^^/M  is  the  mass  ratio  of  the  mineral  contained 

N 

in  the  rock.  It  is  evident  that  Z  a  =1  and  a  are  the  per- 

n  n  n 

tentages  in  Table  V.  Equation  (17)  may  be  applied  at  each  pres~ 
sure  or  stress  level  of  interest,  and  may  also  be  thought  of  as 
representing  a  mass- weighted  averaging  technique. 

2.  Direct  Synthesis  Method 

As  may  be  seen  from  equation  (17),  an  adequate  test  of  th..  ; 
simplified  synthesis  model  requires  a  knowledge  of  the  shock¬ 
loading  characteristics  of  the  rock,  its  mineral  constituents, 
and  the  minerals'  mass  fractions,  a  .  For  many  of  the  materials 
for  which  Hugoniot  data  have  been  published,  mineral  or  petro¬ 
graphic  analyses  sufficient  for  our  purposes  have  not  been  made. 

In  some  cases  it  is  necessary  to  resort  to  the  average  lis.ing 
of  mineral  percentages  (Table  V).  Use  of  equation  (17)  requires 
Hogoniot  data  to  be  in  terms  of  specific  volume  and  pressure. 

Pres sure -volume  Hugon5ot  data  are  shown  in  Table  VI  for  several 
composites  and  in  Table  VII  for  several  minerals.  The  values  of 
a^  used  for  the  Shoal  granite  are  taken  from  reference  26.  These 
values  of  the  mass  fractions  a^  specify  granodiorite  rather  than 
granite. 

A  schematic  representation  of  the  synthesis  of  an  arbitrary 
material  is  provided  in  the  third  figure  of  reference  27.  The 
results  of  the  synthesis  method  using  equation  (17)  and  the  values 
of  the  compressed  volumes  from  Table  VII  are  shown  in  Figure  33. 
The  mass  fractions  a^  used  in  the  synthesis  are  based  on  those 
given  in  references  26  and  28  for  granodiorite  and  basalt,  respec¬ 
tively.  These  values  were  modified  to  include  only  those  minerals 
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Table  VI 


HUGONiOT  DATA  (COMPRESSED  SPECIFIC  VOLUMES)  FOR  IGNEOUS  ROCKS 
(USED  IN  EQUATION  20  FOR  RESULTS  IN  FIGURE  34) 


Preficure 

(kilobars) 

Granite 
(Ref.  10) 

Granodiorite 
(Ref.  12) 

Gabbro 
(Ref.  29) 

Dunite 
(Ref.  29) 

300 

0.246 

0.258 

0.246 

mm 

250 

0.266 

0.276 

0.260 

■n 

200 

0.286 

0.294 

0.275 

0.268 

150 

0.306 

0.313 

0.289 

0.276 

100 

0.327 

0.332 

0.304 

0.286 

50 

0.348 

0.352 

0.318 

0.296 

0 

0.372 

0.372 

0.333 

0.308 

Table  VII 


HUGONIOT  DATA  (COMPRESSED  SPECIFIC  VOLUMES)  FOR  SEVERAL 
MINERAL  CONSTITUENTS  OF  IGNEOUS  ROCKS 
(USED  IN  EQUATIONS  17  AND  20  FOR  RESULTS  IN  FIGURES  33  AND  34) 


Pressure 

(kilobars) 

()uartz 
(Ref. 2, 3) 

Orthoclase* 

Oligoclase* 

labradori te* 

Olivine* 

300 

0.248 

mem 

0.251 

0.232 

tm 

250 

0.268 

wBim 

0.271 

0.254 

mm 

200 

0.288 

0.291 

0.277 

mm 

150 

0.309 

C.314 

0.311 

0.300 

0.268 

100 

0.323 

0.335 

0.332 

0.324 

0.280 

50 

0.342 

0.359 

0.354 

0.346 

0.298 

0 

0.375 

0.390 

0.377 

0..,69 

0.312 

“Volumes  for  these  materials  were  taken  from  the 
early  results  of  experiments  in  this  program 
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nthesized  Hugoniot 

goniot  data  for  gabbro  (Ref.  29) 

gon'ot  data  for  basalt  (Ref.  28) 


Lgure  33(b)  Comparison  of  Synthesized  Hugoniot  with  Experimental  Data  for  Gabbro 
Basalt  (Synthesis  based  on  mineral  Hugoniot  data  of  labradcrite  and 


Synthesis  Hugoniot  data 
Hugoniot  data  for  dunite  (Ref.  29) 
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Figure  33(c)  Comparison  of  Synthesized  Hugoniot  with  Experimental  Data  for  Dunite 
(Synthesis  based  on  mineral  Hugoniot  data  for  olivine) 


listed  in  Table  VII  for  which  Hugoniot  data  exist.  In  addi¬ 
tion,  the  relationship  Za^  =  1.0  was  satisfied.  The  experimentally 
obtained  Hugoniots  for  the  composites  synthesized  are  also  shown 
in  Figure  33.  Two  sets  of  Hugonit t  data  are  shown  for  basalt  and 
gabbro.  One  set  corresponds  to  that  of  reference  28,  the  other 
to  that  of  reference  29. 

3.  Indirect  Synthesis  Method 

Hugoniot  synthesis  may  be  accomplished  in  a  somewhat  indirect 
manner  in  cases  where  only  incomplete  mineral  data  arc  available 
but  partial  composite  data  are  also  available.  To  take  advantage 
of  this  partial  knowledge,  equation  (17)  may  be  generalized  as 
follows ; 

^i (composite)  ~  ^ 

t"  h 

where  V.  is  the  volume  of  the  i*""  composite  and  a.,  is  the  mass 
r  IJ 

f"h  i'h 

fraction  of  the  j  mineral  in  the  i*"  composite.  The  summation 

hh 

of  the  partial  volumes  is  over  the  j  minerals  in  the  i*"  composite. 
Equation  (13)  may  be  applied  at  any  pressure  so  that,  with  further 
generalization. 


Z  a.  .V., 

j  "J 


(19) 


th 

where  k  refers  to  the  k^  pressure  so  that  refers  to  the 
th  th 

volume  of  the  i*"  composite  at  the  k''  pressure.  Equation  (19) 
may  be  expanded  into  matrix  form  for  clarification. 

In  equation  (20),  each  row  represents  a  Hugoniot,  the  com¬ 
posites  in  the  lefthand  V-matrix  and  the  minerals  in  the  v-matrd.x. 
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Composite  Hugoniot  Mass  Fraction  Data  Mineral  Hugoniot 
Data  Data 


Since  some  of  each  are  known  (solid  arrows),  and  some  unknown 
(dashed  arrows),  the  equation  is  rewritten  and  solved  for  the 
un  known  Hugo  n io  t  s . 

This  procedure  was  carried  out  for  the  materials  listed  in 

Table  V.  The  values  of  a„  were  taken  from  Table  V  but  modified 

n 

as  shown.  The  results  of  this  calculation  are  shown  in  Figure  34. 
Since  Hugoniot  data  for  these  materials  do  not  exist,  these  in¬ 
direct  calculations  cannot  serve  as  a  check  or  test  for  the  syn¬ 
thesis  method  but  rather  as  an  illustration  of  its  ultimate 
application. 

4.  Temperature  Calculations 

To  gain  insight  into  the  question  of  the  validity  of  the 
assignment  of  an  average  or  effective  temperature  to  a  mineral 
composite,  temperatures  were  calc\ilated  along  the  mineral  Hugoniots 
and  on  several  unloading  isentropes.  Details  of  this  assignment 
and  possible  modifications  to  the  synthesis  method  are  included 
in  Appendix  B.  The  methods  described  in  reference  18  were  used 
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Normal  pressure,  P  (kilobars) 


Figure  34  Synthesized  Composite  Hugoniot  Curves  Calculated 
from  Equation  24  Using  Data  of  Tables  VI  and  VII 
and  Modified  Values  of  Table  V 


i 

and  will  not  be  detailed  here.  A  computer  program  was  written 
t:o  calculate  the  required  temperatures.  This  calculation  in¬ 
volved  the  solution  of  the  equations  given  in  reference  18  for 
any  two  of  the  four  Hugoniot  variables,  P,  V>  I'g  and  Up. 

The  results  of  die  temperature  calculation  are  given  in 
Table  VIII.  These  results  are  plotted  in  Figure  35.  These  calcu¬ 
lations  include  results  for  granite.  It  is  interesting  to  note 
that  the  granite  results  fall  in  the  mid-range  of  temperatures 
for  the  higher  pressures.  The  various  granites  listed  are  des¬ 
cribed  in  Figure  35.  There  are  two  reasons  for  the  slight  nega¬ 
tive  temperature  change  in  one  granite  sample.  First,  the  computer 
program  utilized  the  rough  rather  than  smoothed  Hugoniot  data.  A 
line  integral  involved  in  the  calculation  is  sensitive  to  abrupt 
or  apparent  changes  in  the  slope  of  the  input  Hugoniot.  Secondly, 
the  method  of  reference  18  is  best  suited  to  Hugoniots  having 
very  low  yield  points,  such  as  the  metals  for  which  it  was  used 
in  reference  6,  or  to  gas-like  materials  having  Hugoniots  with  no 
cusps. 

The  temperature -volume  Hugoniot  data  are  shown  in  Figure  36. 
Unloading  adiabats  were  also  calculated.  These  unloading  paths 
lie  above  the  Hugoniots.  Several  such  adiabats  are  listed  in 
Table  IX.  The  methods  described  in  Appendix  B  may  be  more  signi¬ 
ficant  to  the  unloading  process  than  to  the  loading  process  becaxise 
of  the  longer  times  involved.  The  physical  constants  used  in  ap¬ 
proximate  temperature  calculations  are  listed  in  Table  X. 

C.  Summary  of  Synthesis  Method 

Comparison  of  the  synthesized  Hugoniots  with  those  determined 
by  direct  experiment  (Figure  33)  indicates  that  the  method,  as 
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Normal 


Figure  35  Temperatures  on  the  Hugoniots  of  Granodiorite 
and  its  Major  Mineral  Constituents 

Hugoniot  data  for  granodiorite  ■))  are  taken 
directly  from  reference  12,  while  data  for 
granodiorite  (2)  are  taken  from  Figure  30. 

The  thermodynamic  properties  used  in  the 
calculation  are  listed  in  Table  X. 


Figure  36  Temperature  Volume  Data  for  Granodiorite  and  its  Major  Mineral  Constituents 


applied,  has  promise.  The  results  are  somewhat  inconclusive  ^ 

because  of  the  uncertain  values  for  the  composite  H^ijoniot  data 

as  well  as  for  the  values  of  the  mineral  mass  fractions  a  .  Never- 

^  1 

theless,  these  results  are  perhaps  as  good  as  would  be  expected,  j 

i 

i 

considering  that  the  synthesis  equation  used  (equation  (17))  is  j 

basically  static  in  nature  and  that  the  mineral  Hugoniot  data  I 

are  limited,  incomplete  and,  iii  some  cases,  questionable. 

t 

The  basalt  curve  as  shown  in  Figure  33(b)  was  synthesized 
on  the  basis  of  the  mineral  analysis  given  in  reference  28.  Be-  | 

cause  of  lack  of  Hugoniot  data,  it  was  not  possible  to  base  the  | 

synthesis  on  all  the  materials  given  in  that  analysis.  In  fact,  j 

the  synthesis  was  based  on  only  two  (labradorite  and  olivine)  of  J 

the  five  materials  given.  Furthermore,  the  data  for  this  material  I 

show  a  relatively  high  initial  volvime,  indicating  the  presence  of  I 

voids,  which  are  not  taken  into  account  in  the  analysis.  It  I 

s 

should  be  noted,  however,  that  the  Hugoniot  data  approach  the  I 

synthesized  curve  at  higher  pressures;  at  these  pressures,  the 
voids  are  probably  compressed  out.  The  Hugoniot  data  for  gabbro  _ 

is  also  shown  and  seem  to  agree  much  better  with  tlie  synthesized  J 

curve.  The  closer  agreement  in  the  initial  density  between  these  ! 

data  and  the  synthesized  curve  indicates  the  importance  of  the  \ 

initial  density  in  the  mechanical  behavior.  In  Figure  33(c)  the  J 

comparison  is  made  for  dunite,  the  major  constituent  of  v/hich  is  | 

olivine.  Dunite,  synthesized  on  this  basis  (i.e.,  100  per  cent  | 

olivine),  compares  very  well  with  the  dunite  Hugoniot  data.  f 

While  the  matrix  method  yielding  additional  composite  data  ^ 

(Figure  34)  cannot  in  any  way  be  used  as  a  test,  it  does  serve  f 
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as  an  indication  of  one  method  of  application  of  the  method.  The 
data  presented  in  Figure  34  were  determined  without  compensating 
for  the  difference  in  initial  density  brought  about  by  the  afore¬ 
mentioned  adjustments  in  a^.  The  initial  composite  density  is, 
therefore,  in  error  but  the  order  is  correct,  i.e.,  the  coraposite 
with  the  largest  initial  density  remains  so  in  the  synthesis. 

The  fact  that  the  yield  point  is  the  same  in  all  new  com¬ 
pos  it  si  reflects  the  fact  that  the  granodiorite  curve  with  a  yield 
point  of  about  50  kilobars  was  used.  In  this  indirect  method 
(equation  (20)),  input  mineral  Hugoniot  data  are  incomplete  and 
calculated  (partially)  from  the  composite  input  data. 

li.  is  somewhat  reassuring  that  all  of  the  curves  in  Figure 
34  have  the  same  familiar  shape.  To  some  degree  this  result  is 
due  to  the  use  of  a  limited  amount  of  similar  mineral  data.  When 
a  composite  is  made  up  of  only  two  or  three  minerals,  it  is  nec¬ 
essary  to  know  the  associated  mass  fractions  fairly  accurately  as 
errors  in  a^,  for  small  a^  and  n,  can  lead  to  large  errors  in 
the  resulting  calculated  Hugoniots. 


TEMPERATURES  ALONG  HUGOMIOT  CALCULATED  FROM  SHOCK  AND  FREE  SURFACE  VELOCITIES  SHOWN  FOR  ORTSCCLASE 
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Table  X 


THERMODYNAMIC  PROPERTIES  OF  GEOLOGICAL  MATERIALS  USED  IN 

TEMPERATURE  CALCULATIONS 


Specific 

Heat, 


Thermal 

Expansion 

Coefficient 


V 


Material 


“  Cp 


•k 


(ioules/gm) 


* 


(per  -^C) 


Low  Pressure 


Compress  ibilit>j 

Po** 


(cm^/dyne)  Kcm 


K 


dyne 


-a) 


Oligoclase 

Labradorite 

Orthoclase 

Olivine 


0.796  X  10' 
0.766 
0.643 
0.79 


11.2  X  10 

13.2 
15.0 
23.8 


-6 


1.74  X  10' 
1.50 
2.12 
1.74 


•12 


iGranodio- 
rite (1)&(2) 


0.727 


24.0 


2.12 


± 


6.4  X  10^ 
8.8 
7.1 
13.7 


11.3 


V 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  synthesis  methods,  the  work  leading  to  the  application 
of  these  «nd  tfie  results  of  the  application  were  detailed 

in  the  earlier  sections  of  this  report.  In  this  section,  the 
complete  work  of  this  program  is  summarized  briefly.  Some  of  the 
possible  conclusions,  both  optimistic  and  otherwise,  regarding 
the  feasibility  of  applying  synthesis  methods  are  listed.  The 
discussion  of  the  synthesis  method  is  divided  into  two  categories; 

(1)  that  involving  what  was  or  might  have  bean  done 
on  this  and  other,  related  past  programs, 

(2)  that  involving  what  might  be  lone  in  any  future 
program. 

These  are  best  considered  in  that  order. 

This  program  was  actually  considered  to  be  somewhat  of  a 
feasibility  study  of  the  methods  suggested  in  the  excellent  work 
of  Chabai  (Ref.  30)  but  applied  using  mineral  rather  than  ele¬ 
mental  Hugoniot  data.  The  necessity  for  gathering  only  a  minimal 
amount  of  mineral  Hugoniot  data,  coupled  with  the  uncertainties 
in  thase  data,  especially  those  related  to  the  yield  behavior  of 
these  materials,  required  that  the  synthesis  calculations  lack 
the  luxury  of  great  detail.  Even  a  casual  look  at  the  synthesis 
methods  employed  here  reveals  that  the  main  features  of  the  dy¬ 
namic  loading  process,  as  outlined  in  Appendix  B,  are  entirely 
absent.  Nevertheless  the  results  of  even  these  calculations  are 
very  promising  but  not  conclusive  evidence  of  ultimate  feasibility. 
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It  must  again  be  emphasized  that  much  of  the  composite 
Hugoniot  data  that  presently/  exist  are  not  actually  available 
for  detailed  comparison,  in  this  or  future  synthesis  work.  This 
is  due  to  the  unfortunate  circumstance  that  it  has  not  been  cus¬ 
tomary  to  perform  a  detailed  mineral  analysis  on  geological  samples 
before  testing.  It  is  probably  worth  noting  that  much  of  the 
past  Hugoniot  work  on  composites  was  done  to  determine  the  beha¬ 
vior  of  a  particular  composite  and,  in  addition,  to  do  this  on  a 
"crash  program"  basis.  Since  there  is  every  evidence  that  there 
will  be  continuing  interest  in  Hugoniot  data,  the  need  and  de¬ 
sirability  of  2  uependflble  synthesis  method  seems  apparent  and 
especially  so  for  crash  programs. 

It  is  believed  to  be  worthwhile  to  recommend,  then,  the  steps 
necessary  in  the  development  of  an  ultimate  synthesis  model.  In 
general,  the  development  of  such  a  model  must  be  based  on  the  dy¬ 
namic  characteristics  of  the  physical  problem  rather  than  the 
quasistatic  ideas  considered  here.  Such  an  analysis  is  very 
probably  amenable  uo  computerized  averaging  or  synthesis  tech¬ 
niques.  One  of  the  more  promising  aspects  of  the  problem  is 
that  composites  may  be  expected  to  be  characterized  by  bulk  iso¬ 
tropy  even  when  made  up  of  anisotropic  minerals.  In  fact,  in 
similar  elastic  problems  Lame”^  constants  have  been  developed  for 
mixtures  of  anisotropic  materials  (Ref.  31) .  Such  bulk  isotropy 
might  be  expected  to  greatly  simplify  the  over-all  problem.  In 
addition,  if  the  method  is  to  be  extended  to  lighter,  more  por¬ 
ous  materials,  such  factors  as  saturation  and  porosity  must  be 
included.  Some  consideration  must  also  be  given  to  grain  size 
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in  an  ultimate  theory  and  application.  Remarks  concerning  grain 
size  in  Appendix  B  are  believed  to  be  appropriate.  Even  without 
these  significant  features  of  an  ultimate  model,  useful  predictions 
may  be  expected  over  a  limited  range  of  material  properties.  Thus, 
by  considering  the  quasistatic  problem,  as  was  done  here,  the  use 
of  mineral  weighting  factors  may  be  extended  to  other  composite 
materials . 

Before  the  development  of  the  ultimate  problem  is  undertaken, 
several  more  immediate  steps  are  believed  to  be  necessary.  These 
were  suggested  above  previously  and  are  listed  as  follows: 

(1)  A  better  understanding  of  the  behavior  of  the  feldspars 
and  other  minerals  in  the  neighborhood  of  the  dynamic 
yield  point  is  required.  Low-pressure  experiments, 
probably  in  shock  tubes  with  pressures  below  100  kilo- 
bars,  are  best.  Even  before  such  experiments  are  con¬ 
tinued  the  records  such  as  those  discussed  in  Section 
III  should  be  re-evaluated  on  the  basis  of  a  continu¬ 
ous  loading  process  as  demonstrated  in  Appendix  A. 

This  type  of  loading  would  occur  as  a  result  of  a  com¬ 
pression  fan  or  by  the  elastic  reverberations  discussed 
in  Section  III.  These  elastic  reverberations  should 
be  reconciled  to  the  sharp  transitions  found  by  other 
investigators . 

Unfortunately  neither  time  nor  funding  was  available 
to  carry  out  this  envisioned  (computer)  analysis  of  the 
experimental  records  obtained  in  tliis  program.  In  sev¬ 
eral  cases,  limited  analysis  indicated  that  much  more 
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Hugoniot  information  is  available  from  these  records 
than  nas  presently  evolved.  It  is  hoped  that  this  phase 
may  be  completed  in  the  near  future. 

(2)  The  synthesis  techniques  developed  here  should  then  be 
applied  to  this  re-evaluated  mineral  Hugoniot  data. 

Some  consideration  should  be  given  to  the  development 

or  use  of  thermo  gradients  and  conduction  in  the  unload¬ 
ing  adiabat  calculations  and  an  attempt  to  synthesize 
the  unloading  process  should  be  made. 

(3)  As  a  rather  direct  requirement  in  a  slightly  more  ad¬ 
vanced  program,  Hugoniot  data  should  be  obtained  for 
composites,  fabricated  or  natural  but  crnitai ning  only 
minerals  of  known  properties. 

Many  of  the  details  of  these  steps  and  others  of  a  more  com 
plete  analysis  have  been  detailed  previously  (Ref.  31)  and  will 
not  be  repeated.  These  three  specific  recommendations  are  inclu 
ded  here  because  they  arc  direct  consequences  of  the  work  re¬ 
ported  here. 
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Appendix  A 


INDEFINITE  TIME  RESOLUTION 

The  determination  of  the  Hugoniot  depends  on  measurements 
of  the  velocities  of  one  or  more  shock  waves.  In  cases  where  two 
waves  are  transmitted,  it  is  necessary  to  determine  the  free- 
surface  velocity  change  associated  with  the  arrival  and  velocity 
of  the  second  wave.  Assuming  the  simple  t-x  diagram  illustrated 
in  Figure  7 ''a),  this  velocity  change  occurs  abruptly. 

Several  methods  have  been  used  to  determine  the  arrival  time 
of  the  second  wave  by  continuous  or  quasicontinitous  monitoring  of 
the  free-surface  motion.  Two  continuous  recording  optical  methods 
were  considered  for  use  in  this  investigation;  the  reflecting 
wire  method  of  Wackerly  (Ref.  15)  and  the  inclined  mirror  method 
of  Fowles  (Ref,  14). 

In  an  effort  to  compare  these  methods  an  optical  mock-up  of 
these  experimental  configurations  was  set  up  in  the  Ir^boratoiy . 
Results  of  these  sir  'laticn  experiments  and  of  several  fiel"*  ex¬ 
periments  using  the  inclined  mirror  configuration  indicated  that 
considerable  difficulty  would  be  associated  with  the  exact  loca¬ 
tion  of  the  transition  time,  in  jeitheii’/impthod.  It,  therefore,  be¬ 
came  of  interest  to  determine  the  sensitivity  of  the  experimental 
Hugoniot  determination  to  errors  in  the  location  of  the  transition 
time,  and  to  determine  if  records  with  indefinite  transition  times 
could  be  analyzed.  The  graphical  method  used  to  make  the  latter 
determination  and  the  optical  simulation  experiments  are  discussed. 
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Optical  Simulation  Experiments 


The  simulation  set-up  is  shown  schematically  in  Figure  A-1. 

The  free  surface  motion  of  a  sample  was  simulated  by  a  reflecting 
surface  driven  by  an  electric  motor.  The  change  in  reflectance 
of  the  inclined  mirror,  which  occurs  as  the  free  surface  of  the 
sample  comes  into  contact  with  the  mirror,  is  simulated  by  thu 
reflector's  being  driven  into  an  optical  mask.  The  image  of  the 
reflector  mask  interface  was  "hen  focussed  onto  the  adjustable 
entrance  slit  of  the  recording  camera.  An  Avco  360  deg  rotating 
mirror  camera  with  rotor  driven  by  a  slow- speed  synchronous  motor 
was  used.  The  combination  of  camera  rotor  speed  and  reflector 
drive  speed  was  sufficient  to  result  in  angles  of  about  those  ob¬ 
tained  in  some  early  field  experiments.  In  addition,  the  geometric 
arrangement  between  the  light  source  reflector  and  recording  cam¬ 
era  was  approximately  that  of  the  field  experiments. 

A  series  of  runs  was  made  in  which  the  angles  a  and  p  (Figure 
A-1),  as  well  as  the  light  intensity  and  slit  width,  were  varied. 

A  sample  record  is  shown  in  Figure  A-2.  In  this  run,  the  reflec¬ 
tor  drive  motor  was  switched  to  a  higher  speed  midway  in  the  run. 
The  change  in  speed  is  analogous  to  an  instantaneous  change  in 
velocity  of  the  free  surface  caused  the  arrival  of  a  second 
wave . 

In  each  run  a  wire  stretched  across  the  reflector  a  few 
millimeters  from  the  free  surface  was  used  to  simulate  the  re¬ 
flecting  wire  method.  The  wire  trace  is  shown  below  the  mask- 
reflector  interface  in  Figure  A-2.  As  described  in  reference  15, 
there  are  two  images  of  the  wire  seen  on  the  film  record.  One  is 
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(a)  Change  in  free  surface  velocity 
not  determinable  (see  Figure  A-1) 


(b)  Change  in  tree  surface  velocity 
indicated  by  arrows  in  run  3 
and  run  4 

Figure  A-2  Sample  Records  from  Mock-up  Optical  Experiments 
Showing  Inclined  Mirror  Cutoff,  Top  of  (a)  and 
(b),  and  Reflecting  Wire  and  Image  Separation, 
Bottom  of  (a)  and  (b)  (In  (a)  the  free  surface 
motion  and  geometric  variables  are  such  that 
the  change  in  free  surface  motion  (correspond¬ 
ing  to  arrival  of  second  v;ave)  is  not  determin¬ 
able  by  either  a.ethod;  in  (b)  change  may  be  seen 
in  both  traces) 
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due  to  the  wl;:e  Itself;  the  other  to  its  virtual  image  from  the 
reflector.  The  relative  motion  between  these  images  corresponds 
to  twice  the  free  surface  velocity  of  the  reflector  in  these 
experiments,  or  of  the  sample  in  actual  field  experiments. 

The  records  were  compared  for  these  idealized  cases  of  very 
good  reflection  for  the  wire  method  and  distinct  reflector-mask 
interface  for  the  inclined  mirror  method.  There  was  little  dif¬ 
ference  between  the  two  methods  in  the  difficulty  of  locating  the 
transition  in  free  surface  velocity, and  it  was  decided  to  continue 
the  experimental  program  using  the  inclined  mirror  method  since 
that  method  does  not  recuire  a  reflecting  sample  surface. 

Graphical  Analysis  Method 

A  graphical  method  was  used  to  determine  the  sensitivity  of 
calculated  Hugoniot  data  to  the  accuracy  in  location  of  any  tran¬ 
sition  times  in  free-surface  sample  motion.  It  consisted  of  con¬ 
structing  a  t-x  or  free-surface  motion  diagram  from  an  asstimed 
Hugouiot.  The  -lugoniot  used  corresponded  roughly  to  that  of 
granite  but  was  approximated  by  two  straight-line  segments,  as 
shown  in  Figure  A- 3 (a). 

To  construct  the  t-x  diagram,  the  following  equations  were 
used.  The  Indefiniteness  in  the  location  of  the  transition 
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(a)  Assumed  Hugoniot 


P 


(b)  t-x  diagram  calcu¬ 
lated  from  (a) 


X 


(c)  Approximate  t-x  diagram 
showing  arbitrarily 
chosen  segments 


Figure  A-3  Schematic  Representation  of  Graphical  Method 
Used  to  Test  Sensitivity  of  Data  Analysis 
Method  When  Time  Resolution  Is  Insufficient 
or  for  Rounded  Records 
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time,  i.e.,  the  arrival  of  the  second  shock  wave  at  the  free 
surface,  was  imposed  by  replacing  the  calculated  segmented  free 
surface  motion  by  the  curved  line  (shown  dotted  in  Figure  A-3(b)). 
The  new  free  surface  curve  was  then  approximated  by  n  segments 
of  equal  time  increments.  The  related  free  surface  and  shock 
velocities  were  then  measured  graphically,  and  the  P-e  states  cal¬ 
culated  using  the  equations 


“  "n-1 


PN  ~  ^P.N-1^ 


1  -  £ 


N-1 


1 


i=N  (Ugj  -  Upj^) 
i=l  ^^Si  ■  ^P,i-1^ 


A  comparison  of  the  various  n  line  approximations  is  shown  in 
Figure  A-4.  The  construction  indicates  that  the  calculation  pro¬ 
cess  is  convergent  and  that,  for  indefinite  transitions,  analysis 
could  be  made  on  this  basis. 

There  are  several  approximations  involved  in  the  calculation: 
Perhaps  the  most  serious  is  the  inherent  construction  error  asso¬ 
ciated  with  establishing  the  curved  free  surface  diagram  of 
Figure  A-3(b).  Since  this  construction  error  would  not  be  a 
source  of  error  in  actual  record  analysis,  application  of  this 
method  seems  desirable  for  indefinite  or  rounded  records. 

It  is  interesting  to  test  the  method  in  cases  where  the 
material  properties  are  such  that  the  t-x  diagram  would  be  curved 
with  no  sharp  transitions.  Such  a  curved  Hugoniot  is  shown  in 
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Figure  A- 5.  The  resulting  n  segmented  approximations  were  cal¬ 
culated  as  before  and  are  also  shown  in  Figure  A-5. 

Convergence  to  the  actual  Hugoniot  is  more  rapid  than  before; 
thus,  the  analysis  method  is  apparently  applicable  to  rounded 
records  as  well.  In  the  second  example  it  was  necessary  to  approx¬ 
imate  the  Hugoniot  by  several  segments  in  order  to  construct  the 
t-x  diagram.  In  this  example  16  segments  were  used. 

When  many  lines  are  used,  analysis  becomes  tedious  and  quite 
inaccurate.  If  many  records  are  involved,  the  equations  may  be 
programmed  for  electronic  computation,  thereby  eliminating  con¬ 
siderable  construction  and  iteration  error. 
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Figure  A-5  Results  of  Calculations  for  Material  Having  Rounded  Yield  Zone  Show 
Convergence  of  Calculated  (P-e)  States  to  Assumed  (F-t)  States  with 
Number  ox  Segments  or  Waves  Assumed  in  the  Calculation 
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THEORETICAL  MODEL  STUDIES 


Dynamic  Loading 

The  simplified  analysis  leading  to  the  equation 


N 


V  =  E  -?T*  V  =  L  a  V 
c  M  n  n  n 

n  n 


N 

Z 


(B-1) 


is  based  on  the  assumption  that  the  pressure  on  all  mineral  con¬ 
stituents  is  the  same.  This  is  essentially  a  static  rather  than 
dynamic  view  of  the  process.  Even  if  an  equation  such  as  (B-1) 
is  considered  an  averaging  equation  where  the  a^  are  unspecified 
weighting  terms  to  be  determined  from  experiments,  it  is  necessary 
to  assign  an  effective  pressure  which  characterizes  the  composite 
material.  The  existence  of  such  an  effective  or  average  pressure 
is,  in  fact,  necessary  to  describe  the  Hugoniot  state  of  Che 
material. 

To  clarify  the  differences  between  static  and  dynamic  loading 
processes,  it  is  helpful  to  consider  the  layered  material  of 
Figure  B-1.  Here,  we  are  conside"’ng  the  one-dimensional  case 
where  it  is  assumed  that  the  strain  in  a  direction  perpendicular 
to  the  applied  force  is  zero.  In  Figure  B-l(b),  both  materials 
A  and  b  experience  the  same  pressure.  The  total  strain  is  equal 
to  the  sum  of  the  contributions  of  the  two  materials.  This  is 
completely  analogous  to  the  discussion  leading  to  equation  (B-1). 
The  composite  state  would  be  represented  by  P  and  e,j,,  where  e,j.  is 
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oi 


i=A,  B  ''o 


the  total  strain  of  the  composite  material,  and  V  ./V  is  the 

oi  o 

initial  voltime  ratio  for  the  layered  constituents.  For  the  one¬ 
dimensional  case  considered  here,  the  area  perpendicular  to  the 
direction  of  the  applied  stress  is  constant  so  that  the  volume 
ratio  may  be  replaced  by  the  length  ratio. 


Et 


(one-dimensional) 


"i 


and 


V  .  -  V. 
01  1 


oi 


1. 
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Even  for  this  relatively  simple  geometry,  the  case  of  shock  load¬ 
ing  is  very  much  complicated  by  the  multiple  shock  wave  reflec¬ 
tions.  In  Figure  B-2(a)and(b)  the  shock- loading  case  is  considered. 
Here  it  is  assumed  that  a  projectile  is  being  used  to  impose  the 
initial  shock  loading  conditions  on  material  A.  The  resulting 
loading  conditions  are  transmitted  by  the  shock  and  unloading 
waves  as  they  reverberate  through  the  sample.  The  wave  and  par¬ 
ticle  velocities  are  best  described  on  a  phase  or  t-x  diagram. 

Such  a  diagram  combined  with  a  pressure  particle  velocity  (P,  Up) 
plot  depicting  the  possible  loading  (Hugoniot)  and  unloading  states 
of  materials  A  and  B  is  shown  in  Figure  B-2(a)  and  B-2(b),  respec¬ 
tively.  To  facilitate  the  discussion  and  diagrams,  the  following 
simplifying  assumptions  have  been  made: 
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(1)  The  projectile  material  is  also  material  A. 

(2)  The  (P,Up)  loading  Hugoniot  is  a  straight  line, 
corresponding  to  the  case  where  all  loading  is 
below  the  dynamic  elastic  limit. 

(3)  Unloading  occurs  along  the  reflected  Hugoniot 
(i.e.,  Hugoniot  reflected  through  a  vertical 
line) . 

On  the  t-x  diagram  the  lines  represent  velocities  either  particle 
(broken  lines)  or  wave  (solid  for  loading,  dashed  for  unloading). 

Each  time  is  represented  by  a  horizontal  line,  so  that  the  load¬ 
ing  state  can  be  seen  at  a  particular  time  from  the  numbered 
regions.  These  states  are  shown  schematically  on  the  (P,Up)  dia¬ 
gram.  In  Figure  B-3(a)  to  B-3(g),  the  pressure  pulse  is  shown  at 
the  various  times  corresponding  to  the  horizontal  time  intercepts 
of  Figure  B-2(a).  As  may  be  seen  from  Figure  B-3,  the  pressure 
state  is  neither  uniform  in  space  nor  constant  in  time.  As  in 
the  static  one-dimensional  case,  however,  the  volume  is  repre¬ 
sented  by  the  lengths;  the  corresponding  values  of  1^^^,  l^g,  1^, 

1^  and  Ig  are  shown  on  the  diagram. 

The  first  phase  of  the  loading  process  is  complete  at  the 
time  tg.  At  the  time,  t^,  the  first  shock  loading  wave  completes 
its  passage  through  the  sample  (i.e.,  arrives  at  the  right  edge 
of  material  B) .  At  this  time  the  pressure  (state  3)  is  relieved 
at  the  free  surface.  This  is  represented  in  Figure  B-2(a)  by  a 
release  or  unloading  wave  (leftward-sloping  dashed  line  separating 
states  3  and  4).  After  t^,  the  average  stress  (either  compression 
or  tension)  decreases  toward  zero.  As  may  be  seen  from  Figure  B-2(b) 
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and  takes  on  many  of  the  characteristics  of  a  simple  one  component 
material.  Although  the  final  state  in  this  case  is  a  state  of 
zero  pressure  (State  1,  Figure  B-2(b)),  this  is  not  necessarily 
the  case  but  depends  on  the  boundary  conditions. 

Temperature  Differences 

In  most  static  determinations  of  the  equations  of  state  or 
pres  sure- voliime  states,  the  temperature  is  held  constant.  Since 
the  temperature  along  the  Hugoniot  varies,  it  becomes  another 
parameter  which  should  be  considered  in  Hugoniot  synthesis.  The 
quasistatic  derivation  of  equation  (B-1)  neglected  the  temperature 
variation  along  the  Hugoniot, 

As  seen  in  Section  III,  the  temperatures  of  all  the  minerals 
at  any  given  pressure  are  not  the  same.  If,  as  was  assumed  in  the 
synthesis,  all  constituents  are  at  the  same  pressure,  some  thermal 
conductivity  will  occur.  The  resulting  temperature  change  of  any 
constituent  subsequent  to  its  initial  loading  may  cause  some 
thermal  expansion  or  contraction.  The  variation  in  the  thermal 
properties  of  the  various  constituents  may  cause  an  over-all 
change  in  the  volume  of  the  composite.  These  volvime  changes  have 
been  neglected  but  it  is  of  interest  to  see  how  they  might  be 
considered.  In  addition  it  is  of  interest  to  consider  the  validity 
or  meaning  of  the  assignment  of  a  composite  temperature. 

The  Hugoniots  of  a  three- component  composite  are  considered 
in  Figure  B-4.  An  average  temperature  might  be  calculated  sim¬ 
ply  by  averaging  the  temperatures  of  the  constituents  suitably 
weighted,  by  their  mass  ratios,  as  represented  by  the  dashed  line. 

A  slightly  more  physically  meaningful  calculation  might  involve 
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B-4  Schematic  Representation  of  Changes  in  Volume  of 
Three  Minerals  in  Intimate  Contact  but  Shocked 
to  Different  Temperatures  As  Thermodynamic 
Equilibrium  Is  Approached 

Solid  curve  for  minerals  are  those  immediately 
after  shock  loading;  dashed  curves  represent 
a  later  time. 
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an  energy  lalance  between  the  internal  energy  increase  of  the 
composite  and  its  constituents;  thus 


MAE 

c 


m. 


(B-2) 


where  AE  and  AE.  are  the  internal  energy  increases  of  the  com- 
c  1 

posite  and  i’^  constituent  respectively. 

From  equation  (B-2),  an  effective  or  potential  composite 
temperature  might  be  defined  by 

’’’c  =  ^  “l  r  ■'l  ■  I  ^ 

^  i^^c^  i  c 

where  T.  is  the  effective  composite  temperature,  c^  and  c^  the 

specific  heat  at  constant  pressure  for  the  i  mineral  and  the 

composite,  respectively,  and  is  the  ambient  temperature.  To 

establish  thermal  equilibrium,  the  minerals  having  temperature 

T.  >  T  would  tend  to  cool:  those  with  T.  <  T  would  tend  to  be 
1C  1  e 

heated.  A  corresponding  change  in  the  vo3.umes  could  also  be  ex¬ 
pected. 

AVj.  =  a^ATj.  (B-4) 


where  AV^  =  -  Vjj^(H), 

th 

V  (H)  is  the  initial  shocked  volxome  on  the  Hugoniot  of  the  i 
mineral  at  the  pressure  P,  is  the  thermally  expanded  or  con- 
tracted  volume,  and  is  the  i  thermal  expansion  coefficient. 
The  situation  is  schematically  represented  in  Figure  B-4  v^here 
the  various  minerals'  Hugoniots  have  been  saifted  through  AV^j^. 
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Such  "corrections”  have  not  been  taken  into  account  in  the 
synthesis.  Physically,  such  volume  corrections  I'epresent  an 
idealized  situation  where  sufficient  time  is  available  for  heat 
transfer  to  occur.  Although  no  calculations  have  been  made,  it 
is  unlikely  that  the  thermal  conductivities  of  the  minerals  are 
sufficiently  high  or  the  times  involved  sufficiently  long  to 
justify  these  calculations. 
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